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IMPRINT OF INTERGALACTIC SHOCKS ON THE RADIO SKY 

Uri Keshet', Eli Waxman' and Abraham Loeb^-^ 
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ABSTRACT 

Strong intergalactic shocks are a natural consequence of structure formation in the universe. These shocks 
are expected to deposit large fractions of their thermal energy in relativistic electrons {^^ ~ 0.05 according to 
supernova remnant observations) and magnetic fields (^^ ~ 0.01 according to cluster halo observations). We 
calculate the synchrotron emission from such shocks using an analytical model, calibrated and verified based on 
a hydrodynamical ACDM simulation. The resulting signal composes a large fraction (up to a few 10%) of the 
extragalactic radio background below 500 MHz. The associated angular fluctuations, e.g. STi > 260(£,e£,B /5 x 
10"'*)(z^/100 Mnz)""* K for multipoles 400 < I < 2000, dominate the radio sky for frequencies < 10 GHz and 
angular scales I' ^ < 1° (after a modest removal of discrete sources), provided that ^ 3 x 10"'*. The 
fluctuating signal is most pronounced for ly < 500 MHz, dominating the sky there even for £,e(,B = 5 x 10"^. 
The signal will be easily observable by next generation telescopes such as LOFAR and SKA, and is marginally 
observable with present-day radio telescopes. The signal could also be identified through a cross-correlation 
with tracers of large-scale structure (such as 7-ray emission from intergalactic shocks), possibly even in existing 
< 10 GHz CMB anisotropy maps and high resolution ~ 1 GHz radio surveys. Detection of the signal will 
provide the first identification of intergalactic shocks and of the warm-hot intergalactic medium (believed to 
contain most of the baryons in the low-redshift universe), and gauge the unknown strength of the intergalactic 
magnetic field. We analyze existing observations of the diffuse radio background below 500 MHz, and show 
that they are well fit by a simple, double-disk Galactic model, precluding a direct identification of the diffuse 
extragalactic radio background. Modelling the frequency-dependent anisotropy pattern observed at very low 
(1-10 MHz) frequencies can be used to disentangle the distributions of Galactic cosmic -rays, ionized gas, and 
magnetic fields. Space missions such as the Astronomical Low Frequency Array (ALFA) will thus provide an 
important insight into the structure and composition of our Galaxy. 

Subject headings: galaxies: clusters: general — large-scale structure of universe — shock waves — radio 
continuum: general — Galaxy: structure 



1. INTRODUCTION 

The gravitational formation of structure in the universe 
inevitably produced strong, collisionless shocks in the in- 
tergalactic medium (IGM), owing to the convergence of 
large-scale flows. In these shocks, electrons are ex- 
pected to be Fermi accelerated up to highly relativistic 
(> 10 TeV) energies, limited by inverse-Compton cool- 
ing off cosmic mic rowave background (CMB) photons 
iLoeb & Waxman 200 0). The resulting 7-ray emission is 
expected to trace the large-scale structure of the universe. 
Rich galaxy clusters, characterized by strong, high veloc- 
ity accretion shocks, should be de tect ed by future 7-ray 
missions as bright 7-ray sources (p^^oe b & Waxman 2000; 
iTotani & Kitavama 2000: Waxman & Loeb 2 000) in the form 
of accretion rings with bright spots at their intersections with 
galaxy filaments (Keshet et al. 2003). Recently, a possible as- 
sociation of 7-ray radiation (as measured by the Energetic 
Gamma-Ray Experiment Telescope, EGRET) with the loca- 
tions of Abell clusters was identified at a 3(t confidence level 
(|Scharf & Mukherjee 2002, but see also Reimer et al. 2003). 
The integrated 7-ray background resulting from strong in- 
tergalactic shocks was calculated using hydrodynamical cos- 
mological simulations ( Keshet et al. 2003; Miniati 2002), at 
the level of e^dJ /de < 0.15 keV s"' cm"- sr"^ This signal 
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constitutes < 15% of the widely accepted estimates for the 
flux of the extragalactic 7-ray background (EGRB, see e.g. 
Sreekumaret al. 1998; Strong et al. 2003), although a recent 
analysis (Keshet et al. 2004a) suggests that the EGRB flux is 
lower than previously thought, at least by a factor of 2. The 
7-ray background from weak merger shocks was estimated to 
be a factor of > 10 weaker than the background from inter- 
galactic accretion shocks ( Gabici & Blasi 2OOI. 

In addition to the inverse-Compton emission from inter- 
galactic shocks, synchrotron radiation should also be emit- 
ted by the relativistic electrons, as they gyrate in the shock- 
induced magnetic fields. The resulting radio signature is ex- 
pected to trace the structure of the universe at low redshifts 
(z < 1), and to be dominated by rich, young galaxy clusters. 
Indeed, extended radio emission with no optical counterpart is 
observed in about 10 % of the rich galaxy clusters (radio halos 
and radio relics, see Giovannini et al. 1999), and in more than 
a third of the young, massive clusters (with X-ray luminosity 
Lx > 10*^ erg s"', see Feretti 2003). The radio emission has 
been identified as synchrotron radiation from relativistic elec- 
trons, but there are different models for the origin of the elec- 
trons and the magnetic fields involved (for a recent discussion, 
see Baachi et al. 2003). EnBlin et al. (1998) have used obser- 
vations of nine radio relics to suggest an association between 
these sources and structure formation shocks, focusing on the 
possibility that the relics are revived fossil radio cocoons orig- 
inating from nearby radio galaxies, re-energized by diffusive 
shock acceleration. It is interesting to note that recently, large- 
scale diffuse radio emission was discovered around a filament 
of galaxies, possibly tracing an accretion shock on this scale 
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jBaechi et al. 2002^ . 

Waxman & Loeb (2000) have proposed a simple model, 
which allows one to estimate both the radio and the 7- 
ray signatures of intergalactic shocks, produced by electrons 
accelerated from the inflowing plasma. Their model al- 
lows one to estimate the radio and the 7-ray backgrounds 
and their anisotropy characteristics, as well as the signa- 
ture of individual clusters. The model uses dimensional 
analysis arguments to estimate the properties of the virial- 
ization accretion shock of a halo, as a function of redshift 
z and halo mass M. Halo abundance estimates at differ- 
ent redshifts [such as the Press & Schechter (1974) halo 
mass function], may then be used to calculate various ob- 
servables. The Waxman & Loeb model approximates the 
strong accretion shocks as being spherically symmetric, and 
neglects weak merger shocks. A fraction ~ 0.05 of the 
shock thermal energy is assumed to be carried by relativis- 
tic electrons, based on supernovae remnant (SNR) observa- 
tions (for a discussion, see Keshet et al. 2003 1). Observations 
of ^ . 1 iiG magnetic fields in cluster halos ( Kim et al. 1989t 
iFusco-F emiano et al. 1999; Rephaeli et al. 1999) require that 
a fraction ~ 0.01 of the shock thermal energy be transferred 
into downstream magnetic fields. With these assumptions, 
Waxman & Loeb (2000) found that strong fluctuations in the 
predicted radio signal seriously contaminate CMB anisotropy 
measurements at sub-degree angular scales and frequencies 
below 10 GHz. 

Identification of radio or 7-ray emission from intergalactic 
shocks holds a great promise for advancing current knowl- 
edge of shock formation in the IGM. It should provide the 
first direct evidence for such shocks, revealing the underlying 
large-scale cosmological flows. When combined with 7-ray 
detection, the radio signal will provide a direct measure of 
the unknown magnetic fields in the IGM, possibly shedding 
light on the processes leading to IGM magnetization. Emis- 
sion from large-scale shocks traces the undetected warm-hot 
IGM at temperatures 10^ K < T < 10^ K that is believed to 
con tain most of the b aryons in the low redshift universe (see 
e.g. lDave et al. 

lool- 

Moreover, the signal may be used to 
study non-thermal physical processes in the intergalactic en- 
vironment, such as Fermi acceleration in low density shocks. 

In this paper we explore the radio signature of intergalac- 
tic shocks and examine its observational consequences. Our 
results derive from a generalized version of the model of Wax- 
man & Loeb (2000), adapted for a ACDM universe, modified 
to incorporate non-spherical accretion shocks, and calibrated 
using the global features of a simulated ACDM universe in 
a hydrodynamical cosmological simulation. We perform an 
analysis of the radio sky in order to assess the feasibility of 
observing the predicted signal; in this context, we review var- 
ious known foreground and background radio signals, and an- 
alyze observations at low (v < 500 MHz) frequencies, where 
the diffuse extragalactic background is most pronounced with 
respect to the Galactic foreground. Our results have impli- 
cations for existing high resolution radio telescopes (e.g. the 
Very Large Array ); for next generation ground-based radio 
telescopes such as the LOw Frequency Array (LOFAR^) and 
the Square Kilometer Array (SKA^); and for future space- 
borne telescopes such as the Astronomical Low Frequency 

* see http://www.vla.nrao.edu ?as'tro/guides/vlas| 

' see http://www.lofar.org 

' see http://www.skatel escope.orgl 



Array (ALFA^). 

In ^ we calculate the radio signal from intergalactic 
shocks. We start by generalizing the model of Waxman & 
Loeb (2000, illustrated in Figure [0, and adapting it for a 
ACDM universe. The non-uniformity of the thermal injection 
rate along the shock front is incorporated into the model by 
modifying the modelled shock front area. The free parame- 
ters of the model are calibrated using global features (such 
as the average baryon temperature) of a simulated ACDM 
universe, according to a previously analyzed hydrodynamical 
cosmological simulation ( Sprin gel et al. 200 1 ). The predic- 
tions of the calibrated model are then shown to agree with 
the results of the cosmological simulation, regarding the ra- 
dio emission (se e Figures|2]|3land|5|i and its 7-ray counterpart 
(see Keshet et al. 2003). We estimate the energy fractions 
and S^B using SNR and cluster halo observations, and evaluate 
the uncertainty of these parameters. A qualitative agreement 
is shown to exist between observations of galaxy-cluster radio 
halos, and model predictions (see Figurel^. 

In |3lwe examine the observational consequences of the ra- 
dio signal for present and future radio telescopes. We assess 
the contribution of various foreground and background sig- 
nals to the brightness of the radio sky (see Figure|6j and to the 
angular power spectrum (APS, see Figure^). In particular, we 
examine the synchrotron foreground from our Galaxy, and the 
contamination from discrete radio sources. We also calculate 
the point-source removal threshold required so that the fluctu- 
ations in the intergalactic shock signal dominate the angular 
power spectrum on small angular scales. Possible methods by 
which the signal may be identified using present and future 
observations are discussed. Other extragalactic radio signals, 
namely bremsstrahlung emission from Lya clouds and 21 cm 
tomography, are also reviewed. 

In l^we analyze the diffuse low frequency (v < 500 MHz) 
radio background (LFRB). After highlighting important ob- 
servational features (see Figure |8}, we present a model for 
the Galactic foreground that allows us to (i) try to sepa- 
rate between the Galactic foreground and the extragalactic 
background; (ii) examine if a simple Galactic model can ac- 
count for the observed Galactic foreground; and (iii) demon- 
strate the importance of observations in the frequency range 
1 MHz < < 10 MHz, where absorption in our Galaxy is sig- 
nificant. We show that existing observational data is consis- 
tent with a simple double-disk Galactic model (see Figure|9|l. 
The implied lack of direct evidence for a diffuse extragalactic 
radio background (DERB) is discussed, given earlier models 
for the LFRB and unsubstantiated claims for direct identifica- 
tions of the extragalactic component. Finally, we show how 
future observations at very low (1-10 MHz) frequencies, e.g. 
by the ALFA, may be used to disentangle the distributions of 
Galactic cosmic -rays, magnetic fields and ionized gas. 

Finally, 13 summarizes our results and addresses their po- 
tential implications. We discuss the conditions under which 
emission from intergalactic shocks could be identified, taking 
into account various present and future telescope parameters, 
confusion with competing signals, and possible systematic er- 
rors in our model. Some consequences of a future positive 
detection of the signal are mentioned. 

2. MODEL 

Here we study the extragalactic radio signal expected from 
^ see |http://sgra.jpl.nasa.gov/htinl_dj/ALFA.html| 
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Table 1. Cosmological model and structure formation parameters. 



Parameter 


Meaning 


Value 


h 


Hubble parameter 


0.67 


k 


Curvature 







Matter energy density 


0.30 




Dark matter energy density 


0.26 




Baryon energy density 


0.04 


Ua 


Vacuum energy density 


0.70 


X 


Hydrogen mass fraction 


0.76 


n 


Fluctuation spectrum slope 


1 


0"8 


Spectrum normalization 


0.9 



the strong intergalactic shocks involved in structure forma- 
tion. We begin by reviewing the dime nsion al-analysis based 
model of Waxman & Loeb (2000) in Next, the model 

is generalized, and adapted for a ACDM universe. In ^2.21 
we modify the halo mass function used by the model, and 
incorporate the spectral features of the emitted radiation. 
In ^2.3l we calibrate the unknown dimensionless parameters 
of the model, using a hydrodynamical ACDM simulation 
iSpringel et a l. 20 01). The energy fractions and are 
eva luate d in ^2.41 and their uncertainty is estimated. Finally, 
in ^2.51 we compare some observational features of galaxy- 
cluster radio halos to the predictions of the model. 

We show that after adjusting its parameters to pro- 
duce some global features of the simulated universe, 
the model yields radio and 7-ray signals consistent 
with those extracted independently from the simulation 
(iKeshet et al. 2003; Keshet et al. 2004b, in preparation). The 
calibration scheme incorporates the inhomogeneous thermal 
energy injection rate along a shock, e.g. in the form of 
'hot spots' at the intersections with galaxy filaments. This 
effectively enhances the synchrotron (but not the inverse- 
Compton) signal, because the magnetic energy density is 
higher in these regions. In addition, both synchrotron and 
inverse-Compton power is shifted to smaller angular scales. 
The integrated spectrum is shown to be nearly flat, with ;//^ 
peaking at frequencies near 100 MHz. 

In the following we use a 'concordance' ACDM model of 
Ostriker & Steinhardt (1995) - a flat universe with normal- 
ized vacuum energy density Ha = 0.7, matter energy density 
Qm = 0.3, baryon energy density J7b = 0.04, Hubble param- 
eter h = 0.67, and an initial perturbation spectrum of slope 
n = 1 and normalization = 0.9. The various parameters of 
the cosmological model are summarized in Table[2 Note that 
our cosmological model is slightly different from the ACDM 
model used by Waxman and Loeb (2000). 

2.1. Dimensional Analysis 

Using dimensional analysis arguments, Waxman & Loeb 
(2000) have related the mass M of a virialized halo, its ve- 
locity dispersion cr, its smooth mass accretion rate M through 
strong shocks, its temperature T and the typical radius r^/, of 
its accretion shock, by 



M-- 



V2 a\M,z) 
5 GH(z) 



M(M,Z)=fa, 



'.2.5 X 10-'7„„/z7ofl-'/V«)^yr" 



(1) 



(2) 



T(M,z) = fTk-s'^laHM,z) 

~ 1 .8 X lOVr h%^a-^giaf^^M%^ K , (3) 



and 



r.,(M,z) = /.-^^ 

~ 1 .9frh-f/\g(a)-^^'M\i' Mpc . (4) 

Here H(z) = IQQ h a'^^^ g(a) km s"' Mpc"' is the Hubble pa- 
rameter, with a = ( 1 + z)"' and g(a) = [il„, + H^a^ + ( 1 - f^m - 
f7A)a]i/2, ;i~0.65 nip is the average mass of a particle (in- 
cluding electrons), M14 =M/10''*Mq, and ks is the Boltz- 
mann constant. For consistency with previous work, we use 
the notation /i7() = /z/0.70 and scale the results using 7 ilb/flm 
(~ 0.93 in the present model). The free parameters face, fr 
and /,. are dimensionle ss fa ctors of order unity, which are as- 
sumed (and shown in ^2.3> to be roughly constant in the red- 
shift range and cosmological model of interest, and must be 
calibrated separately. The definitions of the model parameters 
are summarized in T able | 2| along with their values as cali- 
brated in ^23land in ^241 

Note that the above estimates refer to strong shocks only; 
weak shocks that may result from mergers of comparable 
mass objects are ignored since they lead to accelerated elec- 
tron distributions with littl e energy in highly r elativistic elec- 
trons (for a discussion see lKesh^ t ^t al, 2 (^Q^[). For example, 
merger trees have been used ('Gabici & Blasi 20031) to show 
that the 7-ray background from merger shocks is a factor 
of > 10 lower than the background from accretion shocks. 
Moreover, this emission is dominated by the strong shocks 
associated with mergers between halos of mass ratio > 100, 
where the distinction between accretion and merger is vague 
(see e.g. Salvador-sole et al. 1998). Weak merger events may 
only enhance the predicted radio background from intergalac- 
tic shocks, and only in v ery low photon frequencies. We shall 
revisit this issue in ^2.21 when discussing the radio spectrum. 

Collisionless, non-relativistic shocks are known to accel- 
erate electrons to highly relativistic energies. Keshet et al. 
(2003) have shown that SNR observations suggest that the 
energy density of relativistic electrons, accelerated by a struc- 
ture formation shock, constitutes a fraction ~ 5% of the 
thermal energy density behind the shock (up to a factor of 
~ 2, see also ^2.41 for discussion). The maximal energy an 
electron can be accelerated to is limited by cooling, predom- 
inantly through inverse-Compton scattering of background 
CMB photons, yielding a maximal e lectron Lorenz factor 
7„,„, -3x10^ ( Loeb & Waxman 2000). 

As stated above, the model focuses on strong shocks, which 
accelerate electrons to a power-law distribution of index p = 2 
in the differential number of accelerated electrons per elec- 
tron energy (equal energy per logarithmic interval of electron 
energy). The luminosity of a halo of mass M at redshift z due 
to inverse-Compton scattering of CMB photons, may thus be 
estimated as 



i/L'^(M,z) = 



1 



21n7„,fl.i- 



ilb M(M,z) 



^e^kBT(M,z) 



1.1 X 10**^ ( facch — 
y acc J I Q 

X a-''\(afl^] M^f erg s 



,5/3 



(5) 
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Assuming that the energy density of the downstream magnetic 
field constitutes a fractio n ~ 1 % of the downstream ther- 
mal energy density (see ^2.4> . implies that the magnetic field 
strength is given by 



B{M,z) = O.U 



It 

fi 0.01 



1/2 



'70 



1/2- 



consistent with observation of galaxy cluster halos (see 
for discussion). This yields a synchrotron luminosity 

BiM,zf/%Ti 



(6) 
l24l 



vLT{M,z)-. 



Ucmh(z) 



■yLf{M,z) 



1.9 X 10 



39 



face fr 



6 



/2 0.05 0.01 



(7) 



-5/2g(fl)13/3 



M 



.7/3 



ergs 



where Ucmh is the energy density of the CMB. Note that 
whereas the inverse-Compton signal depends on the param- 
eter only logarithmically (through the maximal energy at- 
tained by the relativistic electrons), the synchrotron emission 
scales almost linearly with ^b- 

Given the number density of halos of a given mass at a given 
redshift, dn(z)/dM, one may integrate equation (jS) or equa- 
tion Q to predict the inverse-Compton or the synchrotron 
background from intergalactic shocks, 



dz- 



cdt 



dM 



dn{z) vL^{M,z) 



(8) 



dz J ""' dM 47r(l+z)4 ' 
For example, the Press-Schechter halo mass func- 
tion gives, for our cosmological model (summa- 
rized in Table [Q, an inverse-Compton background 
flux (uP^) = 1 .3 faccfr (6/0-05) keV s"' cm'^ sr"' , 
and a synchrotron background flux {vl'^'") = 3.2 x 
10-i%./r/r'(e./0.05)(eB/0.01) erg s"' cm'^ sr"'. The 
halo mass function, the emitted radiation fields and the 
integration procedure, are all illustrated in Figure^ 

Waxman and Loeb (2000) have calculated the two-point 
correlation function of the radiation emitted by intergalactic 
shocks. The low optical depth of large, hot clusters, which 
dominate the background, enables one to neglect cases where 
more than one cluster lies along the the line of sight, and ap- 
proximate 

5^IA^)={I,{n)I,{y))-{I,f 



dz 



cdt 
dz 



dM 



1 



dn(z) 
dM irrshiMjZ)^ 



vL^{M,z) 
47r(l+z)4 



xP, 



1|2 



i/"^a(z) 



(9) 



r,h{M,z) 

where u and v are unit vectors that satisfy u • v = cost/;. The 
function Pi|2 is the probability that one line of sight passes 
through a halo of radios r.,/,, given that another line of sight 
passes through the same halo, 

Pil2(x)=- f del dyyQ[l-{y+xco?,ef-{x?,\n6f] 



l-±^ 
27r 



X V4--x^+ 2 arcsin(ji:/2) 
}e(2-x). 



+ 2 arctan(;c/ ^A-x^) 



(10) 



where 9 is the heaviside step function. Halo-halo correla- 
tions, neglected in equation (|9}, can only enhance the two- 
point correlation function. The fractional correlation func- 
tion, ^("0) = \/ 5'^Iu{ip)/ {K), is independent on the dynam- 
ical parameters fj, face, £,e and £^b, and depends on the 
geometrical factor /, only through the relation** ^y{ijj',fr) = 
fr^£,v{fr^4'', !)■ The contribution of halos of different mass 
and redshift to 5^1^ is illustrated in Figure [2 



2.2. Adaptation for ACDM , Spectral features 

The Press-Schechter halo mass function is known to dis- 
agree with numerical simulations, predicting less rare, mas- 
sive halos and more abundant, low mass halos. Sheth, Mo, & 
Tormen (2001) have shown that the Press-Schechter approach 
agrees better with cosmological simulations, if the naive crit- 
ical over-density Sdz) indicating collapse (and calculated for 
spherical collapse) is replaced by 



SAz,M) = V^Sc(z)\ l + b 



a\M) 



aSc(z)^ 



(11) 



where (j{M) is the variance of the density field smoothed 
on a mass scale M, and a, b and c are dimensionless 
parameters of order unity. Sheth et al. have shown that 
a modification of this form is naturally obtained if one 
considers ellipsoidal collapse and the effect of shear, al- 
though this requires that a = 1 and does not explain the 
deficit in massive halos, suggesting the importance of ad- 
ditional effects. Numerically, good agreement with ACDM 
simulations is obtained if one chooses a = 0.73, b = 0.34, 
and c = 0.81 (Jenkins et al. 2001; Barkana et al. 2001). 
Such a modification enhances the expected background 
from intergalactic shocks, because it increases the density 
of rare, massive halos, which dominate the extragalac- 
tic signal (Figure ^ illustrates the difference between 
the mass functions). We thus find a 20% increase in 
the predicted inverse Compton background, (i^/jf) = 
1-6/flcf/r (Ce/0.05) keV s"' cm"^ sr"'. The synchrotron sig- 
nal is even more sensitive to the abundance of massive halos, 
and is enhanced by the above modification by 50%, {vPJ") = 
4.9 X IQ-'^faccflf;^ (6/0.05) (6/0.01) erg s"! cm'^ sr"' 
(see Figures n and |2j. In what follows, we use the modified 
halo number density, unless otherwise stated. 

The spectrum of the integrated radiation from a halo ac- 
cretion shock is essentially a broken power-law, where the 
spectral break is introduced by (inverse-Compton dominated) 
cooling, and the break frequency i^briM, z) is related to the age 
(or equivalently, the formation redshift) of the halo. We focus 
on strong shocks which accelerate electrons to a power-law 
energy distribution with an index p = 2. Hence, at frequencies 
below the spectral break the spectrum scales as cx j/'/^, 
reflecting the original distribution of electrons accelerated by 
the shock, whereas above the spectral break, cooling results 
in a flat (j//^ cx i^") spectrum. The spectral break frequency 
roughly corresponds to the minimal energy, at which elec- 
trons manage to significantly cool (say, by 77 = 50% of their 
initial energy) between the cosmic time of their host halo ac- 
cretion shock (given by its redshift z) and the present epoch 

* At low (u < vtr) or high {v > Vmax) frequencies (see i|2.2l . the scaling 
relation is more complicated. 
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Fig. 1 . — Illustrating the model of ^in the redshift — halo mass plane. 
Shown are the halo number density (upper-left panel, color scale: log|Q {rfn(z,M)/rfln(A/) [Mpc"']}) and the fractional halo contributions to three signals: 
inverse-Compton emission vl'^ (upper-right, logjQ/,c [eV s"' cm"^ sr"']), synchrotron emission vl^" (bottom-left, logjQ/,v,, [erg s"' cm"^ sr"']), and syn- 
chrotron squared fractional two-point correlation function ^u(ip)^ (bottom-right, log|o/j2). Each of the three signals is obtained from the corresponding fractional 

halo contribution / by integration, according to F = J dz J d ln(M) f. The contribution of halos in a given part of the M-z phase space may thus be evaluated, 
whereas integration over the entire depicted range of M and z will produce results close to the total signal predicted. 

The results shown coiTespond to the modified halo mass function i Sheth et al. 2001 : Jenkins et al. 2001 ; Barkana et al. 2001 , see i|2.2l . with non-calibrated free 
parameters fr = face = fr = 1. and energy fractions = 0.05 and = 0.01. The dash-dotted curves delineate the relatively small regions, for which fine tuning 
the halo mass function has lowered the halo number density. Spectral features are shown for some photon frequencies for inverse-Compton (contour labels in 
log[Q e/keV) and for synchrotron emission (labels in logjQ u/MHz), as solid lines where i/ coincides with u,„ax{M,z), and as dashed lines where v coincides with 
vi,r(M,z)- The two-point correlation function is shown for an angular separation of i/i = 0.01°; for different tp values, the image is distorted and the angular cutoff 
(bright solid line) shifts (dashed lines, contour labels: i/i/deg). 
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jKeshet et al. 2003h . 



Vhr{M,z)~9h 



70 



1- 



B{M,z) 
0.1 AiG 



{\+zfl^ 



-, -2 



dz 



kHz . 

(12) 

The flat, high frequency part of the spectrum extends up to fre- 
quencies corresponding to the maximal energy of relativistic 
electrons, at which the latter are effectively accelerated by the 
shocks. The maximum Lorenz factor to which electrons are 
accel erated, -fmax ^ 3.3 x 10^(B/0.1 ^lGy^\T/10'^ K)i/2(1 + 
z)~^ (iKeshet et al. 2003h . thus implies maximal synchrotron 
frequencies around 



,-(M,z)~3.6x 10' 



B(M,z) 



0.1 nG 



T(M,z) , 
1^(1+^) Hz, 



(13) 

The spectral features of the emission from intergalactic 
shocks may be shown pictorially in the halo mass-redshift 
plane (see Figure [Q. As seen from equation ( I13t . low mass 
(cool) halos do not contribute to the spectrum at high pho- 
ton frequencies, because they do not accelerate electrons to 
sufficiently high energies. Hence, for a given photon fre- 
quency I' we may draw a curve in the M-z plane (the solid 
curves in Figure [!), such that for halos located along the 
curve the cutoff frequency v,„ax{M,z) coincides with v. Ha- 
los with masses lower than found on this curve do not con- 
tribute to the integrated spectrum at v. In a similar fash- 
ion, a given photon frequency v defines a curve in the M-z 
plane (the dashed curves in Figure [0, for which the spec- 
tral break frequency VbriMjZ) coincides with ly. The emis- 
sion from older halos (with redshift higher than found on 
this curve) is flat at e. The resulting synchrotron spec- 
trum, integrated over all halos, is presented in Figure |2] and 
cuts off below ~ 100 MHz. Weak shocks, neglected in the 
model, accelerate softer electron distributions than discussed 
above, so modelling them results in an integrated spectrum 
slightly softer than depicted in the figure (similar to the es- 
timated inverse- Compton spectrum in the 7-ray band, e.g. 
IKeshet etal 2003; Miniati 2002; G abici & Blasi 200 3). Fig- 
ure|2]thus suggests, that for the total emission from intergalac- 
tic shocks (when accounting also for weak shocks), lyl^, peaks 
at frequencies around 100 MHz. 

2.3. Calibration Using a Cosmological Simulation 
We now turn to a hydrodynamical ACDM simulation 



JSpringel et al. 200 U IKeshet et al. 2003h in order to obtain 
rough estimates of the dimensionless free parameters used in 
the above model. We focus on the epoch < z < 2, which 
is most relevant for radiation from intergalactic shocks. The 
halo parameters fj, face, and fr may be estimated by com- 
paring various quantities calculated in the simulation, to their 
values according to the model. For this purpose we use global 
features of the simulated universe, and not the radiation fields 
resulting from the intergalactic shocks extracted from the sim- 
ulation. This enables us to later test the calibration scheme by 
comparing features of the radiation from intergalactic shocks, 
as extracted from the simulation and as calculated from the 
model. Whereas the halo parameters may be calibrated in 
this fashion, the energy fractions and can not be sim- 
ilarly evaluated from a cosmological simulation, but require 
independent observational data (see ^2.4> . When evaluating 
features of the model, we have neglected low mass halos with 
temperatures below r„„„ = 10"* K, where collapse is strongly 



suppressed. Our calibration scheme is insensitive to the exact 
value of T,„i„ . 

The temperature parameter fj may be estimated using the 
temperature statistics of the baryonic component of the uni- 
verse. We examine the mass averaged temperature, well fit 
in the simulation by (T(z))m — 4 x 10^ e""-^-' K for redshifts 
< z < 2. At low redshifts, z < 1, the agreement between 
the simulation and the model is best if fr ~ 0.5. However, 
the redshift dependence of {T{z)) m differs between the model 
and the simulation. For example, agreement between the 
two on the average temperature of the present-day universe, 
{T(z = 0))m, requires fj ~ 0.45; the best fit at the epoch 
< z < 2 is for fr = 0.52; and at the extreme, agreement on 
(r(z = 2)) M requires fr ~ 0.68. These examples imply that fr 
has a weak redshift dependence, and that we should consider 
the range of/7 ~ 0.45-0.55. 

In order to assess the accretion rate parameter face, we ex- 
amine the fraction of mass that has been processed by strong 
shocks between redshift zo and the present epoch, fprodz < 
Zo). According to the cosmological simulation, fprodz < 1) — 
18%, reproduced by the model if /„„ — 0.12. Here too, the 
redshift dependence of the accretion rate differs between the 
simulation and the model, implying a weak redshift depen- 
dence of face- For example, fpmdz < 2) ~ 41% corresponds 
to face — 0.17, whereas the accretion rate at very low redshifts 
implies that face — 0.08. These findings suggest that face lies 
in the range of 0.08 -0.17 for the relevant epoch. 

Using the halo mass function modified for a ACDM uni- 
verse (see ^2.21 and the best fit values found for the above 
calibrated parameters, the model yields an inverse-Compton 
background ulf ~ 0.1(^e/0.05) keV s~' cm"^ sr"', in good 
agreement with the results of the simulation. Note that the 




Fig. 2. — Spectrum of the synchrotron background from intergalactic 
shocks. When using the ACDM-modified halo mass function, the (non- 
calibrated) model predicts a signal (dashed line) with ~ 60% more energy 
flux than found when using the Press-Schechter mass function (dotted line). 

The synchrotron brightness scales as Ii, oc faccfjfr^, so the signal corre- 
sponding to the calibrated free-parameters of Table El ("heavy solid line) is a 
factor of ~ 3 stronger than the signal obtained with non-calibrated parame- 
ters, /r = face = /r = 1. The spect rum calculated from a ACDM simulation 
(dash dotted line.l Keshet et al. 2004h in preparation) considers only photons 
above the spectral break i/i,r{M, z) of each halo. Although the simulated spec- 
trum contains ^ 60% the energy flux of the corresponding spectrum of the 
cahbrated model (when neglecting emission below the spectral break, thin 
solid line), the two spectra are similar at the relevant (< 10 GHz) frequen- 
cies. 
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result is independent of the parameters and ^b- Agreement 
to better than ^ 30% should be regarded as somewhat of a nu- 
merical coincidence, considering the uncertainties and the fact 
that lyl'^ scales linearly with the product frfacc- Nonetheless, 
the good fit suggests that the calibration of the parameters /j- 
and face is sensible. 

Finally, the shock radius parameter fr can be deduced from 
the simulation by studying the morphology of cluster accre- 
tion shocks. A shock ring of diameter 5-10 Mpc has been 
identified around a simulated cluster of mass M(r < 5 Mpc) ~ 
1O'^M0 (Ke shet et al. 2003). Such a range of shock radii is 
obtained by the model, if fr ~ 0.6- 1 .2. Clearly, more work is 
required in order to obtain a better understanding of simulated 
cluster accretion shocks, to better calibrate fr, and possibly 
identify its redshift dependence. 

An important aspect of the model is the sensitivity of the 
predicted synchrotron luminosity of a halo, but not its inverse- 
Compton lumino sity, to halo asymmetry (as well as to ^b, 
as mentioned in ^2.1> . Such asymmetries introduce spatial 
fluctuations in the thermal energy injection rate through the 
halo accretion shock. Whereas the inverse-Compton emis- 
sion from a particular region along the shock front scales lin- 
early with the downstream thermal energy, the synchrotron 
emission scales as the square of the thermal energy, vUJ" cx 
MT^. This implies, that fluctuations along the shock front 
will increase the overall synchrotron luminosity of the halo, 
leaving its inverse-Compton luminosity intact. Accounting 
for fluctuations is crucial, because cosmological simulations 
(e.g.^iniati et al. 2001) find highly asymmetric shocks, with 
very strong fluctuations along the shock fronts. In pa rticular, 
the simulation we have studied (iKeshet et al. 20031) reveals 
"bright spots" in the thermal energy injection rate at the in- 
tersections of the accretion shock of a ~ lO'^M© cluster with 
large galaxy filaments, channelling large amounts of gas into 
the cluster region. The thermal injection rate (and the result- 
ing 7-ray brightness) of these regions is more than an order of 
magnitud e higher than the typical brightness along the shock 
front (see lKeshet et al. 2(X)3t Figures 9 and 10). 

The large fluctuations in thermal injection rate necessitate 
the introduction of a geometrical correction factor into the 
simple model presented in ^2.11 We note the localized na- 
ture of the hot regions along the shock front, and the large 
ratio ~ 10 between the temperature of these regions and the 
temperature of the dimmer, more extended regions. These 
features justify approximating the geometrical correction by 
attributing all the emission from a halo accretion shock to an 
effective smaller region, of scale 7sh — '"ja/C- This implies 
replacing the dimensionless parameter fr with a different pa- 
rameter, fr = fr/C- Note that such a replacement will en- 
hance the synchrotron luminosity of the halo by a factor <^^, 
while leaving its inverse-Compton luminosity unchanged. As 
an order of magnitude estimate, /,• ~ 0.6- 1.2 roughly yields 
fr = fr/C — 0.05-0.20, with a middle-of-the-road estimate of 
fr = 0.1. Note that the geometrical correction alters the two- 
point correlation function of both synchrotron and inverse- 
Compton emission, shifting power to smaller angular scales. 

The model parameters calibrated above are summarized in 
Table |2l We present the resulting synchrotron radio back- 
ground in Figure|2] and the synchrotron two-point correlation 
function in Figure |3] It is encouraging to note that after cal- 
ibrating the free parameters with essentially global features 
of the cosmological simulation, the model predictions agree 
well with results extracted from the simulation, regarding the 



Table 2. Parameters of the intergalactic shock model. 



Param. 


Quantity parameterized 


Value (Range) 


h 


Halo downstream temperature " 


0.50 (0.45-0.55) 


face 


Halo mass accretion rate " 


0.12(0.08-0.17) 


fr 


Halo shock radius " 


0.9 (0.6-1.2) 


fr 


Halo energy injection scale " 


0.1 (0.05-0.2) 




Energy fraction in relativistic electrons 


0.05 (0.02-0.10) 




Energy fraction in magnetic fields 


0.01 (0.005-0.04) 



"Defined in i|2.1l and calibrated in i|2.3l using a hydrodynamical cosmologi- 
cal simulation. 

''Energy fractions out of the shock thermal energy, evaluated in il2.4l accord- 
ing to observations of SNRs and galaxy cluster halos. 



inverse-Compton 7-ray background (IKeshet et al. 2003h . the 
synchrotron radio background, and the synchrotron two-point 
correlation function ( Keshet etal. 2004b). A more accurate 
calibration of the parameters, including an evaluation of their 
redshift dependence, may be obtained from a detailed analysis 
of the various clusters identified in cosmological simulations. 




V [deg] 

Fig. 3. — Fractional two-point correlation function of synchrotron 

background from intergalactic shocks, for v = 100 MHz. With the modified 
halo mass function, the (non-calibrated) model predicts a signal (dashed fine) 
slightly stronger than when using the Press-Schechter mass function (dotted 
line). The result of the calibrated model (dash-dotted line), scaling accord- 
ing to £^v(,^,i';fr) = fr^$,v(fr^ip,yA) for vtr <S <S Umax, has more power 
on small scales. These model predictions are well fit by the functional form 
^j^(ip, v) = A(i/)exp {-[i/)/i/>o(i^)]''*''*}, e.g. a fit is shown (triangles) for the 
calibrated model signal, with A = 25, V'o = 9' and ^ = 1 .5. The ACDM simula- 
tion signal (circles with dashed line to guide the eye, error bars represent sta- 
tistical errors introduced by averaging over random pairs of lines of sight, see 
Keshet et al. 2004h) is similar to the calibrated model result, and yields even 
more small scale power. The above results do not include emission below the 
break frequency vi,r. The results of the calibrated model for the entire spec- 
trum have more large-scale power, due to the contribution of young, nearby 
halos, and are presented for frequencies 1 MHz, 100 MHz and 10 GHz (solid 
lines, bottom to top). 



2.4. Energy Conversion Efficiency 

Intergalactic shock waves are expected to accelerate elec- 
trons to highly relativistic energies, and to strongly amplify 
magnetic fields. The average fractions of shock thermal en- 
ergy deposited in relativistic electrons and in magnetic 
fields (^fi) are both important parameters of our model, each 
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bearing linearly upon the predicted synchrotron signal, and 
thus deserve a special discussion here. 

Electron acceleration 

Collisionless, non-relativistic shock waves are known to 
Fermi accelerate a power-law energy distribution of relativis- 
tic particles. This phenomenon has been observed in astro- 
physical shock waves on various scales, such as in shocks 
forming when the supersonic solar wind collides with plan- 
etary magnetospheres, in shocks surrounding SNRs in the in- 
terstellar medium, and probably also in shocks in many of the 
most active extra galactic sources, quasars a nd radio galax- 
ies (|Drurv 1983* 'Blandford & Eichler 19871) . The electi'on 
power-law distributions extend up to ~ 100 TeV energies in 
SNRs (Tanimo ri et al. 199 8). where shock velocities are of 
order v ^ IQ^ km s~', similar to intergalactic shock velocities. 

Although no existing model credib ly calculates the ac celer- 
ation efficiency, a simple argument JKeshet et al. 20031 ) sug- 
gests that we may evaluate S,e in strong intergalactic shocks 
using the estimated acceleration efficiency of other strong as- 
trophysical shock waves. Consider an ideally strong, non- 
relativistic shock wave, such that the shock Mach number 
T ^ 1 and the thermal energy of the upstream plasma is neg- 
ligible with respect to the shock energy. The physics of such 
a shock is essentially determined by three parameters, namely 
the shock velocity v, the upstream plasma number density n„, 
and the upstream magnetic field strength B„ (in principle, the 
result may also depend on the detailed structure of the up- 
stream magnetic field). The upstream density may be elimi- 
nated from the problem altogether by measuring time in units 
of v~p^ , where Vp is the plasma frequency. The upstream mag- 
netic field strength, parameterized by the cyclotron frequency 
Vc, can not be additionally scaled out of the problem. How- 
ever, comparing Vc^i (where subscript / denotes a property of 
the ions) to the growth rate of electromagnetic instabilities in 
the shocked plasma, z/,-,,, = Vp/iv/c, indicates that their ratio in 
strong shocks satisfies {vd/vins)^ = (B^/87r)/(nj,mpV^/2) ^ 1. 
We thus assume that there is a well behaved limit when this 
ratio approaches zero, implying that the upstream magnetic 
field has little effect on the characteristics of strong shocks. 
With this assumption, we expect to find much similarity be- 
tween sufficiently strong shocks in different environments, 
provided that their shock velocities are comparable, regard- 
less of the plasma density and the strength or structure of the 
upstream magnetic field. The little effect of upstream mag- 
netic fields on strong shocks is supported by recent obser- 
vations of SNR shocks, which suggest that the shocks pro- 
duce strongly fluctuating, near equipartition (with respect to 
the thermal energy of the non-relativistic electrons, hereafter) 
magnetic fields, much stronger than the magnetic field far up- 
stream (Bamba et al. 2003; Berezhko et al. 2003, and discus- 
sion of magnetic field amplification, below). 

The preceding discussion indicates that the best analogy to 
strong intergalactic shocks may be found in the strong SNR 
shocks, drawing upon the similarity between the velocities of 
the two families of shocks, both of order v ~ lO"' km s"' . The 
fraction of shock energy deposited in relativistic electrons by 
SNR shocks was estimated by several authors, with the most 
reliable estimates found for remnants w ith 7-ray detection 
such as SNR1006 (^Tani mori et al. 19981) . Dyer et al. (2001) 
have modelled the multi-frequency emission from SNR 1006, 
finding an acceleration efficiency of = 5.3% (correspond- 
ing to a fraction Ty^ ~ 1 .4% o ut of the total supernova explo- 
sion energy). Other authors jMastichiadis & de Jaeer 19961 



lAh aronian & Atovan 19991) have estimated rj^ ~ l%-2% in 
SNR1006, which corresponds to values in the range of 
3.8% -7.6%. Recent Chandra observations, resolving thin 
sheets in the NE shock front of SNR 1006, suggest electron ac- 
celeration more localized and more efficient than previously 
thought, with near equipartition energy (with respect to the 
thermal e nergy of the non-rel ativistic electrons) in relativistic 
electrons (* Bamba et al. 2003h . The extrapolated value of 
depends on the uncertain magnetic field strength; the lowest 
plausible energy fraction in relativistic electrons corresponds 
to an equipartition magnetic field, for which ~ 4% - 11 % 
for different sheet regions. Ellison, Slane andGaensler (2001) 
have modelled the emission from SNR G347.3 -0.5, estimat- 
ing the energy of relativistic electrons to constitute at least 
1 .2%, and more likely 2.5%, of the shock kinetic energy flux, 
corresponding to2.5%<^f~5%. 

The above estimates of the acceleration efficiency of SNR 
shocks, suggest that for strong intergalactic shocks ~ 0.05, 
with an uncertainty factor of ~ 2. An independent, less reli- 
able method for estimating i^^, relies on the ratio between the 
energies of cosmic -ray electrons and cosmic-ray ions in the 
interstellar medium, sugge sting that ~ l%-3% (for a dis- 
cussion, see lKeshet et al. 2003 ). However, since the relation 
between this ratio in the interstellar medium and immediately 
behind intergalactic (or SNR) shocks is unknown, this esti- 
mate is highly uncertain. We stress that our estimate of 
is based only on observations of SNRs, without provoking 
any elaborate model for the acceleration mechanism of the 
electrons by the shocks. The high efficiency of electron ac- 
celeration deduced, suggests that a substantial fraction of the 
shock thermal energy is transferred into relativistic ions, indi- 
cating that a non-linear theory (e.g. Ber ezhko & Ell ison 19991 
and references therein) is required in order to account for the 
shock structure and the particle acceleration process. Non- 
linear theories for diffusive acceleration of particles by shock 
wav es are at the present stage incomplete (for a recent review, 
see (Malkov & Drurv 20011) . In particular, such models pre- 
dict deviations of the accelerated particle distribution from a 
pure power-law. Our model, which is based on observations 
of SNR shocks, does not incorporate such deviations. 

Magnetic field amplification 

Magnetic fields in galaxy cluster halos have been esti- 
mated using several techniques, based on diffuse synchrotron 
emission from the cluster, preferably combined with inverse- 
Compton detection, based on Faraday rotation of background 
or embedded polarized radio sources, and based on the obser- 
vation of cold fr onts in cluster X-ray images (for reviews, see 
Kron berg 1994t iHenrisen 1998; Carilli & Tavlor 2002). The 
synchrotron emission from a cluster measures the volume- 
averaged magnetic field weighted by the relativistic electron 
distribution, and is thus most relevant for our model. Studies 
of the resulting radio signal suggest volume averaged mag- 
netic field strengths of the order of a few 0.1 /iG, close to the 
values inferred for a minimal energy configuration with equal 
energy in magnetic fields and in relativistic electrons. Fara- 
day rotation and cold front studies suggest stronger magnetic 
fields of several /iG, but are sensitive to different measures 
of the magnetic field. Faraday rotation studies, for example, 
suggest magnetic fields of several /iG with a high filling factor 
up to ^ 0.5 Mpc from the cluster center (Clarke et al. 2001). 
However, such studies are sensitive to the magnetic field 
weighted by the thermal electron distribution, suffer from an 
uncertain level of depolarization internal to the source, and 
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their analysis depends on the assumed configuration of the 
magnetic field. 

In the Coma super-cluster, as the most studied exam- 
ple, magnetic fi eld measuremen ts based on the cluster- 
size radio halo (Kim et al. 1989*), excess emission in the 
X-ray (Renhaeli et al. 1999; Fusco-Femian o et al. 1999) and 
in the extreme UV (Hwan2 1997, but see also Bowyer et 
al. 1999) ban ds, and lack of detection in the 7-ray band 
(jSreekumar et alJ99^, all suggest a volume-averaged mag- 
netic field of a few 0. 1 /iG, close to the value of ^ 0.4 /iG ob- 
tained by assuming equal energy in magnetic fields and in rel- 
ativistic electrons ( Giovannini et al. 1993). Studies of Fara- 
day rotation of background polarized sources, on the other 
hand, suggest magnetic fields of several /iG, ent angled on 
< 1 kpc scales ( Kimetal. 199 0-. Feretti et al. 1995h . 

A fraction — 0.01 of a shock thermal energy transferred 
to magnetic energy, implies magnetic field strengths that are 
an order of magnitude lower than their equipartition value, 
and are close to their value in a configuration with equal en- 
ergy in magnetic fields and in relativistic electrons (assum- 
ing ~ 0.05). In cluster halos, our calibrated model re- 
produces, for this choice of ^b, a volume-averaged magnetic 
field strength B ~ 0.1 ^iG (for M ~ IO'^^Mq), although the 
'bright spots' could contain magnetic field strengths as high as 
1 /iG. The lowest observational estimates for magnetic field 
strengths in such clusters are > 0.05 /iG, suggesting that 
could be higher than our chosen value by a factor of a few, but 
is unlikely to be smaller than it by more than a factor of ~ 2. 

It is important to note that recent studies suggest that SNRs 
contain strong, near equipartition magnetic fields. For exam- 
ple, Pannuti et al. (2003) claim that modelling the broad 
band (radio to 7-ray) observations of SNR G347.3-0.5 re- 
quires a magnetic field strength B ~ 150;!;|o" ^'^'^ implies 
highly localized X-ray emission. The existence of strong, 
near equipartition magnetic fields behind SNR shocks has 
recently been confirmed by the highly localized nature of 
hard X-ray emission in the resolved sheets com posing the 
shock front of SNR1006 (e.g. lBerezhko et al. 2003i) . The nar- 
row extent of emission observed upstream of these sheets 
implies a small upstream diffusion constant, correspond- 
ing to upstream magnetic fields of strength > 10 /iG, with 
strong fluct uations {SB/B ~ 1) on small, d ^ 10'^ cm scales 
JBamba et al. 2003) . Such magnetic field fluctuations are far 
stronger, and have much smaller scales, than found in the sur- 
rounding ISM, and so must be induced by the shock. These 
conclusions strongly suggest that the assumption, that the 
physics of a strong shock depends only weakly upon the far 
upstream magnetic field, is valid, justifying the analogy be- 
tween strong intergalactic shocks and strong SNR shocks of 
comparable velocities. 

2.5. Comparison with Radio Halo Observations 

Radio halos are observed in ^ 35% of the young, mas- 
sive gakx^_clu^ters_(wiAX-^^ Lx > 10*^ erg s"\ 
see iGiovan nini et al. 1999t iFeretti 20031) . and their radio lu- 
minosity is known to be correlated w ith the cluster tem - 
perature (Liang etal. 2000) and mass dGovoni et al. 

lool- 

Figure |3 shows the specific luminosities of such halos for 
V = \A GHz, plotted against the temperature ranges of their 
host clusters. A simple power-law fit to the data, Li,{T,v = 
1.4 GHz) cx T^, gives = 3.55, which is very similar to the 
power-law index = 7/2 predicted by our model, although 
the (calibrated) model yields a specific luminosity that is 
lower than observed by a factor of 8 (equivalently, the tem- 



peratures are discrepant by a factor of ^ 1.8). A similarly 
good qualitative agreement between the observations and the 
model is obtained when replacing the cluster temperatures 
with t heir estimated mas ses, giving L^{M,i/ = 1.4 GHz) cx 
M^^ dGovoni et al. 2 007). similar to the model prediction 
U™(M) (X Nfl^ . However, the estimated mass of a cluster 
is somewhat less certain than the measured temperature, and 
is sensitive to the definition of the cluster boundary. 
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Fig. 4. — Specific luminosities of galaxy-cluster radio halos for v — 
1.4 GHz, plotted against the temperat ure ranges of their host c lusters (error 
bars with cluster names, adopted from, Govoni et al. 200UlFeretti 200 3). The 
best power-law fit to the data (solid line) scales as (T, 1/ = 1 .4 GHz) oc 
T'^^. The specific luminosity according to the model scales similarly, 
L^(r) oc r'''~, but the normalization (of the calibrated model, dashed line) is 
different; the dash-dotted line presents the result of the model where the halo 
temperature was re-scaled by a factor of 1.8. 

Some comments regarding Figure0]are in place. First note, 
that the agreement between model and observations regard- 
ing the value of is insensitive to the free parameters of the 
model. The specific synchrotron luminosity of the model, 
when written as a function of halo temperature, scales linearly 

with the combination faccfr^fj'^^^.eS.B- The strong depen- 
dence of this result upon the model parameters, and the low 
redshifts iz < 0.3) of all the halos depicted in Figure|3 suggest 
that the discrepancy between model and observations may be 
eliminated by better calibrating the parameters, and modelling 
their redshift dependence (note that this will increase the radio 
signal we predict). Second, note that the low fraction of clus- 
ters with observed radio halos could result from a combination 
of two effects: (i) not all clusters have been significantly ac- 
creting new mass at the cosmic time when they are observed; 
and (ii) the surface brightness of the halos observed is low, 
close to the instrumental sensitivity, suggesting that more ra- 
dio halos will be observed with future telescopes. A varying 
accretion rate will tend to enhance the luminosity of the halos 
that are observed, possibly accounting for some of the dis- 
crepancy seen in Figure |4] Finally, note that Miniati et al. 
(2001) have found, using a cosmological simulation, that = 
2.6-2.8 (for emission from electrons accelerated at shocks; 
for emission from secondary electrons produced by p-p colli- 
sions of cosmic ray ions, they found </> = 4.1-4.2). However, 
their simulation differs substantially from our model and from 
the simulation of Keshet et al. (2003, 2004b), regarding their 
cosmological model ( SCDM), electr on acceleration efficiency 
(^c < 0.5%, see also iMiniati 2002^ ■ and magnetic field nor- 
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malization ((B^)'''^ ~ 3 /iG for a Coma-like cluster). 

3. OBSERVATIONAL CONSEQUENCES 

In this section we examine the capability of present-day 
and next generation radio telescopes operated from the ground 
(the LOFAR and the SKA) and from space (the ALFA) to de- 
tect the radio emission from intergalactic shocks. The LOFAR 
is an interferometric imaging telescope for the 10-240 MHz 
band, planned to begin initial operation on 2006. The SKA 
is a square kilometer interferometric array for the 150 MHz- 
20 GHz frequency range, planned to become operational dur- 
ing the next decade. The ALFA mission is a proposed space 
interferometer composed of 16 satellites, for very low fre- 
quencies in the range 30 kHz - 30 MHz. The main features 
of these telescopes are summarized in Table|3] 

This section is organized as follows. In ^3.11 we present 
the notations used throughout the section. Var ious foreground 
and background signals are discussed in ^3.21 in an attempt to 
assess the optimal conditions for detecting the emission from 
intergalactic shocks, with as little contaminatio n as possible. 
We summarize the results of this section in ^3.31 and dis- 
cuss scenarios by which the signal could be identified. A de- 
tailed analysis of low-frequency sky-brightness observations 
and their interpretation will be given in ^ Implications of the 
model uncertainties on the feasibility of detecting the signal, 
will be discussed in ^ 

3.1. Notations 

As mentioned in ^ the two-point correlation function of 
specific intensity fluctuations at a given angular separation ip, 
is defined as 

= {iA»)W))-{i.f , (14) 

where u and v are unit vectors that satisfy u • v = cos ip, and 
Slij(u) = /i/(u)- (/,/). The specific intensity is related to the 
brightness temperature Ti, (which approximately equals the 
thermodynamic temperature for the sky brightness and the 
frequency range of interest) by = Iv^ksTh/c^. The follow- 
ing discussion of intensity fluctuations is thus equally appli- 
cable for temperature fluctuations, up to a multiplicative con- 
stant. 

It is often advantageous to study the Angular Power Spec- 
trum (hereafter APS) of various signals, in order to obtain a 
direct estimate of their importance at various angular scales. 
The angular power spectrum C/ is defined through the relation 

, CO 

S%iij)=—Y,i2l+l)C,(,y)P,(costP), (15) 
1=1 

where P;(x) is the Legendre polynomial of degree /. The 
power at a given multipole I is often expressed using its log- 
arithmic contribution to the intensity variance, Sli = [1(21 + 
l)C//47r]'''^ (or 5Ti for the variance of the brightness tem- 
perature), where multipole / roughly corresponds to angular 
scales ~ 180°//. Using the orthogonality of the Legendre 
polynomials, we may invert equation il5\ to find the APS, 

Ci{v) = 2-K I S^I^iiP)Pi(cosij)smiPd^lj . (16) 
Jq 

The APS of synchrotron emission from intergalactic shocks 
may be calculated from the two-point correlation function 



(see Figure 13 using equation ( I16t . The resulting APS, pre- 
sented in Figure|5j peaks at multipoles I ~ 400-4000 for the 
relevant frequency range, corresponding to angles 9 ~ 3'-30'. 
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Fig. 5. — Angular power spectrum of the synclirotron background from 
intergalactic shocks. The logarithmic contribution to the variance u5Ii = 
!^[/(2/ + l)Q/47r]'/- is shown according to the calibrated model (see i|2.3l 
and Table|2) for frequencies 1 MHz, 100 MHz, an d 10 GHz (dashed lines, 
bottom to top). The result of the ACDM simulation fKeshet et al. 2004b) for 
u = 100 MHz (band, with thickness corresponding to the statistical error), 
neglecting emission below the spectral break, exhibits more power than the 
corresponding result of the calibrated model (when neglecting emission be- 
low the spectral break, dash dotted line), in particular on small, arcminute 
angular scales. 

In order to illustrate the meaning of the APS, we expand the 
fluctuations in the specific intensity by the spherical harmon- 
ics Yi,„(a), 

OO / 

6IAu) = 5m a,,„{iy)Y,,„(a) . (17) 

l=i m=-l 

A radiation field drawn from an isotropic distribution satisfies 
{ai„W))=0 (18) 

and 

Qiv) , (19) 

where () represents averaging over the ensemble. The Legen- 
dre polynomial addition theorem justifies the identification of 
Ci as the APS in equation \\9\ . In practice, for signals which 
are not necessarily i sotropic, such as the Galac tic foreground, 
it is traditional (e.g. lTegmark & Efstathiou 19961) to define 

Ci{v)=:^Y.^\aai^t) ■ (20) 

m=-l 

It is often advantageous to analyze extragalactic signals by 
focusing on a small region P in the sky, where Galactic 
foreground is minimal. When studying angular correlations 
within P, the contribution of multipoles that correspond to an- 
gular scales larger than P is diminished. Hence, for angular 
scales 9 much smaller than P, the correlation function may be 
approximated using equation jl5L by introducing a low mul- 
tipole cutoff at 1{P). For example, the correlation function in a 
sky patch P of typical angular size 15° (~ 0.4% of the sky) in- 
volves only multipoles I > l(P) ~ 10. At the other extreme, a 
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Table 3. Parameters of next generation radio telescopes. 



Parameter 


LOFAR 


SKA 


ALFA 


Frequency range 


10 MHz - 90 MHz 
110 MHz - 240 MHz 


150 MHz - 20 GHz 


30 kHz - 30 MHz 


Angular resolution 


< 10" neai- 15 MHz 

< 1" near 150 MHz 


< 0.1" for 1.4 GHz 


1°.7 near 100 kHz 
1' for 10 MHz 


Continuum surface 
brightness sensitivity 


10 mK on 5' -10' in 
high frequency band 


1 Kon 0.1" 





high multipole cutoff is introduced by the angular resolution 
of the experimental device. Assuming a Gaussian response 
function with standard deviation Oc, this is incorporated into 
equation M51 as an exponential multipole cutoff. With these 
modifications, we conclude that 

'5'UV')=^I](2/+l)Q(i^)^Kcos^)expf-^j . (21) 

i=i(P) ^ ^ 

Finally, the sensitivity of a telescope is often estimated us- 
ing its RMS noise level for a given beam size a = nip'^, de- 
noted Sirms = [<J^/(c)] For well behaved signals, where the 
two-point correlation function S^Iiip) is a monotonically de- 
creasing function of -0, one finds 6^I(ip) < S^Imsi'^)- Hence, 
the RMS noise level sensitivity imposes an upper limit to the 
two-point correlation function sensitivity, and may be com- 
pared directly with the logarithmic contribution to the vari- 
ance at the corresponding multipole, Sli. 

3.2. Competing Signals 

In the following, we examine various foreground and back- 
ground signals in the radio band, for the relevant frequency 
range (1 MHz "^v <\Q GHz). We begin with Galactic syn- 
chrotron foreground and with the integrated signal from dis- 
crete radio sources, which contaminate the expected fluctu- 
ating signal from intergalactic shocks on large and on small 
angular scales, respectively. We then discuss other extragalac- 
tic radio signals of interest, namely bremsstrahlung from Lya 
clouds and 21 cm tomography. For each signal discussed, we 
present its contribution to the sky brightness in Figure |6l and 
its logarithmic contribution to the variance on a 0° .5 angular 
scale in Figure 

Galactic Synchrotron Emission 

At low frequencies v <\ GHz, the brightness of the sky 
is dominated by Galactic synchrotron emission, produced by 
cosmic -ray electrons gyrating in the magnetic fields of the in- 
terstellar medium (ISM). At low Galactic latitudes, there is 
also a substantial contribution of free-fre e emissio n from low- 
latitude H II regions (e.g Bacciaalupi et al. 200 Ih . In the fre- 
quency range of interest (v < 10 GHz), the two-point correla- 
tion function is dominated by Galactic synchrotron emission, 
because of its significant power on large angular scales. The 
observed brightness of the radio sky on these scales, domi- 
nated by Galactic foreground, is discussed extensively in ^ 
Here we focus on the level of contamination by Galactic syn- 
chrotron emission in high Galactic latitudes and on small an- 
gular scales. 

With large uncertainties regarding the spatial distribution 
of Galactic cosmic rays and magnetic fields, the best esti- 



mates of Galactic synchrotron emission are obtained by ex- 
trapolating direct measurements carried out at frequencies and 
angular scales where synchrotron emission is reliably mea- 
sured. We thus make the key assumption, often used in 
the literature (in par ticular in CMB anisotropy studies, e.g. 
'Tegm ark et al. 2000l) . that the multipole dependence of the 
Galactic synchrotron APS varies little with frequency, such 
that Ci{v) ~ f{l)g{v). The frequency dependence g{v) of 
Galactic synchrotron emission is discussed in ^ The APS 
multipole dependence /(/) may be extracted from high reso- 
lution maps at ^ GHz frequencies down to arc-minute scales 
(multipoles I ~ 6000, Tucci et al. 2002), and may in principle 
be extrapolated into the entire frequency regime and multipole 
range of interest. 

However, the angular power spectra extracted from radio 
maps depend on a number of factors, most importantly the 
Galactic latitude range examined and the efficiency at which 
discrete radio sources are removed from the map. Unfortu- 
nately, high resolution surveys are available mostly for low 
latitudes, which are less promising as potential search sites for 
an extragalactic signal because they exhibit a stronger Galac- 
tic foreground, in particular on small angular scales. Nonethe- 
less, improved high and medium latitude data have recently 
made it possible to extend our knowledge of the high latitude 
APS up to multipoles I ~ 800. 

The high latitude (60.5° < b < 84.5°) low resolu- 
tion (FWHM 0.6° - 2.3°) maps of Brouw & Spoelstra 
(1976) for 5 frequencies between 408 MHz and 1411 MHz, 
have been used to estimate the APS up to / ~ 70 (e.g. 
BruscoH et al. 2002, table 6). Analys is of the Bonn 408 MHz 
full-sky survey ( Haslam et al. 19821) suggests a power-law 
APS of the form Ci{l > 1) oc I'P with /? ~ 2.5-3.0 down 
to the survey resolution limit 0.85°, corresponding to multi- 
poles / ~ 200 (TegrnaA & Efttathioii 1996)- With a charac- 
terist ic coherence sc a le ^ 5°- 10° for synchr otron emission 
(e.g. i Spoelstra 19841 iBandav & Wolfendale 1991.) . Tegmark 
& Efstathiou (1996) argue that Q cx (/ + 5)"^ describes the 
APS rather well. Thus, although the contribution of Galac- 
tic synchrotron emission to the two-point correlation function 
is significant, it is dominated by large angular scales around 
5° - 10° (low multipoles I < 40) and introduces little con- 
tamination on small angular scales (see Figure 0. The APS 
exhibits strong fluctuations across the sky, suggesting that 
'quiet' regions may be identified and selected as preferable 
search sites for an extragalactic signal. 

Analysis of the 2.3 GHz Rhodes map (Uonas et al. 1998ft . 
with FWHM resolution 20', reveals strong variations of 
the APS with Galactic latitude, and a large contribution 
of discrete rad io sources which tend to flatten the APS 
( iGiardino etal. 20 01). At high latitudes (\b\ > 20°), Giardino 
et al. (2001) find, after removing discrete sources using me- 
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dian filtering, that 

Q ~ AH-f^ , (22) 

with A = 0.3 ± 0.2 K, and /3 = 2.92 ± 0.07, valid up to / ~ 100. 
Although the quoted uncertainties are large, comparison with 
the low resolution studies mentioned above indicates that the 
normalization A is probably no larger than 0.3 K. Similarly 
steep APS, with (3 varying in the range 2.60-3.35 for dif- 
ferent medium latitude {\b\ < 20°) regions, were found from 
a 1.4 GHz survey carried out with the Effelsberg 100 m tele- 
scope ( Uvaniker et al. 1999), a fter remov ing > 5a unresolved 
discrete sources (Baccigaluni et al. 2001 ). The high angular 
resolution of the Effelsberg survey, 9c ~ 9.35', thus suggests 
that the steep power-law found by Giardino et al. holds up 
to multipoles I ~ 800. Note that after removing the brightest 
intensity peaks associated with low latitude H II regions, the 
APS highly resembles the APS of the polarized intensity com- 
ponent (.Baccigalupi et al. 2001) . as expected for synchrotron 
emission. 

Discrete Radio Sources 

Discrete radio sources, mostly radio galaxies, active galac- 
tic nuclei (AGNs) and normal galaxies, make an impor- 
tant contribution to the extragalactic radio background. 
Catalogues of radio sources at several frequencies have 
been used to estimate their contribution to the radio sky 
iSimon 19771 ^illis et al. 1977). This procedure is limited 
by the uncertain contribution of faint, unidentified sources 
(mostly normal g alaxies), which domi nate the background 
(see for example Windhorst et al. 1993) and must be mod- 
elled. Other studies have estimated the contribution of dis- 
crete sources to the radio background, utilizing the well- 
known correlation between the infra-red and radio flux den- 
sities of individual galax ies (Protheroe & Biermann 1996; 
iHaarsma & Partridg e 19981) . However, the results of such 
studies are sensitive to the unknown redshift evolution of the 
sources. In addition, the radio-infra-red correlation holds only 
in cases where the radio emission is associated with star for- 
mation. Hence, such studies must be supplemented by inde- 
pendent estimates of sources not associated with star forma- 
tion, such as a cata logue-based estimate of the emission from 
AGNs (Rv leT968h . 

The studies outlined above suggest that the contribution of 
discrete sources to the radio sky is roughly an order of mag- 
nitude lower than the observed sky brightness. This corre- 
sponds to 25% -65% of the total extragalactic background, 
as estimated b y spectral modellin g of low frequency radio 
observations ( C lark etal. 1970- Cane 19791 s ee ®, although 
some evolutionary models (e.g. lProtheroe & Biermann 1996) 
may even account for the entire strong extragalactic signal 
calculated by Clark et al. (1970). The spectral index of 
the integrated emission from discrete sources has been esti- 
mated to lie in the range s = 0.7-0.8 for frequencies in the 
range 100 MHz -few GHz JSimon 19771 iLawson et al. 19871 
IHaarsma & Partr idge 1998V although some dependence upon 
scale may be expected because the source distribution is es- 
sentially bimodal (ITegmark et al. 2 000). At very low (j^ < 
10 MHz) frequencies, the background from discrete sources 
turns around, probably because of synchrotron self-absorption 
JSimon 

It is easier to estimate the contribution of discrete sources 
to the APS, because the latter is dominated by the bright, well 
studied sources. Discrete radio sources have approximately 
a Poisson distribution in the sky, because projection through 



their wide redshift distribution effectiv ely dim inishes their 
correlations (e.g. Il^gmark & EfstathiouT^S^- Hence, as 
long as the angular scales concerned are much larger than the 
angular extent of the sources, the two-point correlation func- 
tion vanishes and the APS is flat, Q cx f (white noise), such 
that the logarithmic contribution to t he variance is 6li ocl. O ne 
may estimate the APS as (e.g. Teg mark & Efstathiou 19961) 

C, = S^-g^dS, (23) 

where S is the source flux and dN/ dS is the differential num- 
ber density of discrete sources in the sky. The limited sensi- 
tivity S,„i„ of the observations limits our knowledge of dN/ dS 
to sources with flux S > S„,m, but since the number density of 
faint sources is not too steep, the possible error introduced is 
small. The upper limit of the integral, essentially determined 
by the brightest sources, can be lowered in order to reduce the 
noise, by modelling and removing the brightest sources from 
the analyzed map. However, this quickly becomes laborious 
as the number of sources increases, and introduces inevitable 
errors associated with source removal uncertainties. 

Following Tegmark & Efstathiou (1996), we use source 
counts produced by the 1 .4 GHz VLA FIRST all-sky sur- 
vey. A small fraction (^ lO"'*) of the discrete sources in 
the FIRST catalogue have angular scales larger than an arc- 
minute (White et al. 1997), so the use of equation i23\ is le- 
gitimate when dealing with multipoles / < lO'*. The number 
density of radio sources (after performing a r esolution correc- 
tion for extended sources. Wh ite et al. 1997^ is well fit in the 
flux range S = I mJy - 1 Jy by 

1^ . 1.5 X 10« f ^V'-Vl.^V Jy-' sr-' . 
dS \lm]yj \ lOOmJyy ^ 

(24) 

For faint sources in the range S = 10/iJy- 1 mJy, the num- 
ber count steepens to dN/dS oc S~^ with 7 = 2.2 ± 0.2 
jWindhorst et al. 19931) . This changes the integrated bright- 
ness considerabl y, bu t has only a negligible effect on the APS. 
Using equation 923\ . we find that for a modest cut at S„,a.x = 
100 mJy, corresponding to removal of the ~ 7 x 10"* brightest 
sources all-sky, Qiiy) ~ 1.7 x 10"'^(i//1.4 GHz)""^ Us - 
ing the 6C survey earned out at 151 MHz (iHales et al. 1988h . 
we find that the above cut is equivalent to removal of sources 
brighter than ^ 3 Jy at 150 MHz. As shown in Figure0 such 
a cut places the emission from discrete sources at the same 
level as the Galactic foreground for angular scales ^ 0°.5. At 
smaller angular scales, contamination from discrete sources 
becomes increasingly worse, requiring the removal of more 
sources. 

Bremsstrahlung Emission from Lya Clouds 

Loeb (1996) had calculated the thermal bremsstrahlung 
emission from Lya clouds. These clouds are most likely in 
photo-ionization equilibrium with the UV background, per- 
mitting an estimate of the bremsstrahlung signal if the UV 
background intensity, 7"' (z), is known. In the following, we 
update the predicted bremsstrahlung signal according to re- 
cent UV observations. 

Based on the neutral hydrogen density estimated from Lya 
absorption lines in the spectra of background QSOs, the emis- 
sion from Lya clouds with column densities below the Lyman 
limit (A^ < 10'^ cm"^) at a redshift z < 5, is given by: 

~ 10"^^-^=^° '*(y22> erg s"' cm-2 sr"' Hz'^ . (25) 
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Fig. 6. — Brightness of the radio sky. Ground-based (circles and error bars for the Galactic polar regions) and space-based (triangles for IMP-6 spacecraft 
minimum and maximum) observations are dominated by Galactic synchrotron emission (dash dotted line shows the modelled Galactic foreground toward the 
polar regions, see ^ in low frequencies, and by the CMB (double-dashed line) at high frequencies. Various estimates of the extragalactic component (legend, 
see text) are typically an order of magnitude lower than the Galactic foreground. The background from intergalactic shocks according to the caUbrated model 
(of ^^and Tabled solid line) is roughly of the same magnitude, and scales according to oc faccfjfi^ ■ The signal is also shown according to the ACDM 
simulation (accounting only for emission above the spectral break, dotted line, see Keshet et al. 2004b). The frequency ranges of the LOFAR and the SKA are 
shown (double arrows and corresponding shaded regions), as well as constant brightness temperature contours (dotted lines, labelled on the top and on the right 
axes). 
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Here (722) was defined as the UV intensity above the Lyman 
limit (wavelengths A < 912A) weighted over its redshift evo- 
lution 



{J22) = 



_ /o^22(z)(l+Zr<'-^** 



(26) 



where J22 = 7"'(z)/(10"^^ erg s"' cm"^ sr"' Hz"'), and the fit 
parameter 7 = 1 .5 ± 0.4 is determined by th e redshift de pen- 
dence of the H I cloud surface density (see iLoeb 1996 1. and 
the references therein). The signal is flat for frequencies 
V <^ 10''* Hz, persisting throughout the frequency range of 
interest. 

The UV background intensity may be estimated using the 
proximity effect, whereby one compares the effects of photo- 
ionization by the UV background and by a nearby QSO on a 
Lya cloud. Scott et al. (2002) have thus estimated the UV 
background in redshifts z < 5 as follows: 



Jiiiz) - 
implying that 

= 10" 



6S+'' ** 

1 0+3.8 
^■^-0.22 

7 0+3-4 
. '-^-4.4 



for z < 1 ; 
fori <z< 1.7; 
forl.7<z<5. 



22 ±1 



erg s 



-1 -2 

cm sr 



Hz" 



(27) 



(28) 



This result neglects the contribution of optically thick clouds 
at high iz > 5) redshifts, and possible additional ionization 
fields. Hence, the actual bremsstrahlung flux may be some- 
what higher than calculated in equation ( I28> . 

Loeb (1996) had also evaluated the amplitude of intensity 
fluctuations introduced by the finite number of H I clouds 
along any given line of sight. Assuming a random cloud dis- 
tribution in redshift and in column density, one finds 



AJ^ = 10 



-23±0.4 



(J'i 



22)'/^ ergs"' cm"^ sr"' Hz"' 



where we have defined 



,2. JoJiiizfil+zT'^'-^^dz 



{J22) ^ 



(29) 



(30) 



/(,(l+z)"(''-7)c/2 

Using equation ( I27> . we thus find AJ^, = 
]^Q-23±i gj-g g-i gjjj-2 jj2"'. Such fluctuations are 
expected on ^ 7" scales (corresponding to I ~ 10^), although 
weaker fluctuations introduced by clouds associated with 
large-scale structure may appear even on ^ 10' scales. This 
fluctuation level should be regarded as a lower limit, because 
in addition to the assumptions leading to equation ( I28> . the 
non-randomness of the cloud distribution further enhances 
the fluctuation signal. Nevertheless, without a unique spectral 
or temporal signature, confusion with the other signals 
discussed in this section, in particular discrete sources, is 
likely to preclude detection of bremsstrahlung from Lya 
clouds throughout the frequency range of interest. 

IGM21 cm Tomography 



Madau, Meiksin & Rees (1997, see also lTozzi etal. 20001) 
have proposed that the 21 cm (1 .4 GHz) spin flip transition of 
atomic hydrogen could be used to probe the IGM at the epoch 
prior to reionization. Spatial inhomogeneities in the IGM may 
be observed today as redshifted emission or absorption fluc- 
tuations against the CMB, in situations where early sources of 
radiation decoupled the IGM spin temperature from the CMB 



temperature. In principle, the combined angular and spectral 
signal, stronger than the CMB fluctuations by two orders of 
magnitude, may be used to trace the 'cosmic-web' structure 
of the early universe, in both space and cosmic time. 

However, as pointed out by Di Matteo et al. (2002), at the 
relevant frequency range 50-200 MHz, contamination by dis- 
crete radio sources imposes a serious contamination even for 
the maximal signal amplitude AT ~ 10 mK expected. Nev- 
ertheless, a spectral ^10 mK feature, caused by the fast rise 
of the Lya background as the first UV sources reionized the 
IGM, could possibly be detected behind the spectrally contin- 
uous foreground (Zaldarriaga, Furlanetto, & Hernquist 2003; 
Morales & Hewitt 2003; Loeb & Zaldarriaga 2003; Gnedin 
& Shaver 2003, and references therein). Emission from in- 
tergalactic shocks provides an additional important source of 
confusion for 21 cm tomography, but its smooth spectrum 
may not prevent detection of the sharp 21cm spectral features. 

3.3. Summary and Conclusions 

We conclude that the design of next generation radio tele- 
scopes such as the LOFAR, the SKA and the ALFA, is more 
than sufficient for detection of the angular fluctuations intro- 
duced by intergalactic shocks, as calcul ated in ^and con - 
firmed by a cosmological simulation ( Kesh et et al. 2004bl) . 
Identification of the signal is limited by confusion with Galac- 
tic foreground and with discrete radio sources. Foreground 
fluctuations in the synchrotron emission of our Galaxy con- 
strain a clear detection of the signal to sub-degree scales 
and to high Galactic latitudes. Confusion with discrete radio 
sources requires that the brightest sources be modelled and re- 
moved. With a feasible point source cut (100 mJy at 1 .4 GHz, 
or 3 Jy at 150 MHz), the signal dominates over the competing 
signals at angular scales 10' < < 1° (by roughly an order 
of magnitude according to the ACDM simulation), whereas 
detection on arcminute scales will require a more ambitious 
point source removal. The spectrum of emission from inter- 
galactic shocks indicates that the signal is most pronounced 
at 100 MHz frequencies, planned to be covered by both the 
LOFAR and the SKA. We find that other extragalactic signals 
in the radio band, namely bremsstrahlung from Lya clouds 
and angular (but perhaps not spectral) 21 cm tomography, will 
be confusion limited by intergalactic shocks as well as by ra- 
dio point sources. 

Interestingly, detection of the fluctuating signal from inter- 
galactic shocks is just possible with present high resolution 
radio telescopes, such as the VLA (at the maximal sensitivity 
configuration, see FigureQ. The calculated signal should be 
detectable at high latitudes, a few 100 MHz frequencies and 
sub-degree scales. As noted by Waxman & Loeb (2000), the 
signal may have akeady been detected by CMB anisotropy 
studies, at ^ 10 GHz frequencies and 9 <\° scales. In par- 
ticular, the signal may have been detected in arcminute scales 
by low frequency CMB anisotropy experiments such as the 
Australian Telescope Compact Array (ATCA), Ryle, and the 
VLA. It is possible, that so me of the noise removed from the 
corresponding maps (e.g. Subrahm anvan et al. 200 0) is asso- 
ciated with emission from intergalactic shocks, and is thus 
correlated with tracers of large-scale structure. 

The high resolution of present ^ GHz high resolution sur- 
veys raises the possibility that their maps already include 
the signal at an identifiable level. Such surveys are mostly 
available at low latitudes (say, \b\ < 8°), and thus exhibit 
stronger Galactic foreground with flatter angular power spec- 
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Fig. 7. — Logai'ithmic contribution SIi of vaiious signals to tlie variance at multipole / = 400 (corresponding to 6 ~ 0°.5). The bacJcground from in- 
tergalactic shocks is shown according to the calibrated model (of i|2.3l and Table l2l solid line), and scales roughly according to SIi(l', £,e,£,B,facc,fT , fr) = 
£,e(,Bfaccfjfr^ ^Ilifrl', 1,1,1,1,1)- The signal according to the ACDM simulation (diamonds with statistical error bars for / = 400, 10^ and lO", from bottom to 
top, accounting only for emission above the spectral break, see Keshet et al. 2004b) scales linearly with ^e^B- The high-latitude Galactic foreground (dash-dotted 
line) and the integrated emission from discrete radio sources (dotted line, assuming removal of sources above San = 100 mJy at 1.4 GHz) are presented. At fre- 
quencies above ~ 10 GHz, CMB fluctuations take over (dashed line extrapolated from WMAP measurements, eiTor bars present detections by the Cambridge Ra- 
dio Telescope on 0° .2—0° .5 scales, and by the Owens Valley Radio Observatory on 0° . 1 —0° .6 scales; for references see Hinshaw et al. 2003^.White et al . 1999). 
Bremsstrahlung emission from Ly a clouds (low horizontal shaded region) and the maximal 2 1 cm tomography signal (dashed contour) are seen to be much weaker 
than the above competing signals. Continuum surface brightness sensitivities of the LOFAR (star) and the SKA (square), and the side-lobe confusion limit of the 
VLA (triangles, for 10 minutes in configuration D, see http://www.vla.nrao.edu/astro/guides/vlas i are shown (labels denote beam widths). Constant brightness 
temperature contours are shown as dotted lines, labelled on the top and the right axes. Hmission from intergalactic shocks is seen to dominate the sky at this 
angular scale for frequencies u < 500 MHz, provided that ^e^B ^ 10^*. According to the simulation, the emission dominates the sky at this frequency range for 
e<0°.5, evenif^(,gs = 5x 10"^ 
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tra than found in high latitudes. This implies significant fore- 
ground contamination, especially on small angular scales. 
Some low resolution studies have found an APS power in- 
dex ~ 2.4 extend ing up to / ~ 900 (Giardino et al. 200^ 
iTpgmarket al. 2000), others finding /3 ~ 1.7 extending up to 
/ ~ 6000 (Tucci et al. 2002). In addition, detection of the sig- 
nal becomes increasingly difficult for frequencies away from 
the optimal frequency, ~ 100 MHz. Nevertheless, large lon- 
gitudinal fl uctuations have been id entified in the low latitude 
APS (e.g. Baccigalu pi et al. 200 Ih . suggesting that the signal 
may surface in the 'quiet regions'. Analysis of such regions, 
including careful removal of discrete sources, may yield the 
desired signal, identifiable at ~ 10' scales by cross-correlating 
the maps with known tracers of large scale structure. 

4. ANALYSIS OF THE LOW FREQUENCY RADIO SKY 

Next, we analyze the diffuse low frequency {v < 500 MHz) 
radio background (LFRB), dominated by Galactic syn- 
chrotron emission. A simple model for the Galactic emission 
is presented in order to (i) try to separate between the Galactic 
foreground and the extragalactic background; (ii) examine if a 
simple Galactic model can account for the observed Galactic 
foreground; and (iii) demonstrate the importance of observa- 
tions in the frequency range 1 MHz < ^ 10 MHz, where 
free-free absorption in our Galaxy is significant. 

We find that the observations are well-fit by a simple 
double-disk Galactic model, given the existing observational 
uncertainties. This implies that the uncertainties and low reso- 
lution of present low frequency observations preclude a direct 
identification of the diffuse extragalactic radio background 
(DERB). Upper limit on the DERB flux at frequencies < 30 
MHz may be imposed, at a level roughly an order of magni- 
tude lower than the Galactic foreground. Point sources and 
intergalactic shocks akeady account for a substantial fraction 
of this upper limit, suggesting that it is not far from the true 
extragalactic signal. 

We point out the lack of models accounting for the com- 
bined spectral and angular data at very low (1 MHz ^ < 
10 MHz) frequencies, where absorption in our Galaxy is non- 
negligible. We discuss the frequency-dependent anisotropy 
pattern observed in this frequency range, reflecting different 
emission and absorption along different lines of sight. We 
demonstrate, when constructing the Galactic model, how the 
combined spectral and angular data may be used to disentan- 
gle the distributions of Galactic cosmic-rays, magnetic fields 
and ionized gas, leading to an elaborate three-dimensional 
model of the Galaxy. Future high resolution observations at 
very low frequencies, produced for example by the ALFA 
space mission, could thus provide valuable information re- 
garding the structure and composition of the Galaxy. 

We begin by reviewing various spectral and angular features 
of LFRB observations, in ^4.11 A simpl e mo del for Galactic 
emission and absorption is presented in ^4.21 and shown to fit 
the combined spectral and angular data rather well. Previous 
LFRB models, which in general account for either the spec- 
tral or the angular aspects of the observations, are reviewed 
in ^.31 and claims for direct identification of the DERB are 
discussed. 

4.1. Observations of the Low Frequency Radio Background 

(LFRB) 

The LFRB has by now been measured by a vari- 
ety of ground-based and space-borne instruments. 
Detailed maps of the low frequency radio sky have 



been produced by ground-based experiments and pub- 
Ushed for frequencie s 1.6 MHz - 16.5 MHz_|Ellis 198i 
lEllis & Mendillo 1987^. 22 MHz ('Ro geretal. 1999iK 
34.5 MHz ~(|D warakanath & Shankar 199^1 38 MHz 

(Milo2radov-Turin& Smith 1973), 45 MH z 

(Alvar ez et al. 1997). and 408 MHz JHaslametal. 1982h . 
The radio background in frequencies below 5 MHz is 
inaccessible to most ground-based experiments, because of 
strong absorption by the ionosphere. Notable exceptions are 
the radio telescopes operated in Tasmania, taking advantage 
of the unusually favorable ionospheric conditions prevailing 
in the region (Ellis 1965) to mea sure the background i n 
frequencies as low as 1.6 MHz (lEUis & Mendillo 19871) . 
At even lower frequencies, one must resort to space tele- 
scopes such as the Radio Astronomy Explorer 1 (RAE-1) 
satellite (Cla rk et al. 1970; Alexande r et al. 1 970), the 
radio astronomy ex periment on board the IMP-6 space- 
craft (j Brown 1 973|). th e RAE-2 lunar orbiting satellite 

ovaco & Bro wn 19781) and the WAVES exper iment on 
board the WIND spacecraft ('Manning & Dulk 20011) . Our 
picture of the very low frequency radio sky will be revolu- 
tionized by future high resolution space telescopes, such as 
the ALFA mission, planned to operate at 30 kHz - 30 MHz 
frequencies (see Table|3j- 

Generally, the ground-based experiments achieve better an- 
gular resolution, ranging from sub-degree scales at high fre- 
quencies (e.g. 51' at 408 MHz) to > 10° at low (- 2 MHz) 
frequencies, whereas space-borne telescopes suffer from very 
poor angular resolution, of order ^ 100°. On the other hand, 
ground-based experiments are generally limited not only to 
higher frequencies, but also to smaller fields of view, requir- 
ing a combination of several long-term surveys in order to 
produce an a ll-sky map, such as t he Bonn all-sky radio contin- 
uum survey dHaslam et al. 1982h . Figure|8]illustrates various 
low frequency observations, carried out by both ground-based 
and space-borne telescopes. Ground-based observations were 
extracted from published sky maps, and are quoted for several 
interesting regions on the sky for which multi-frequency data 
exist, such as the Galactic center (GC) and the south Galactic 
pole region (SPR). The data are summarized in Table|3 

Spectra 

The low frequency radio spectrum is often approximated 
as a broken power-law. The specific intensity is found to 
peak at frequencies in the range Vp^ak — 2 - 20 MHz, where 
Vpeak varies across the sky (Figure |8|l. In the frequency 
range 20-200 MHz, the specific intensity I^, measured to- 
wards a given line of sight decreases with increasing fre- 
quency, and is well fitted by a power-law cx v~\ with 
spectral indices in the range s ~ 0.2-0.7 varying across the 
sky. For example, s = 0.55 ± 0.03 is found toward the north 
Galactic pole, s = 0.65 ±0.15 toward the south Galactic pole, 
s = 0.38 ± 0.02 toward the Galactic anti-center (see Figure 
m iBridle 19671 ICane 197*^ and the references therein) and 
s ^ 0.25 towards the Galactic center At low frequencies 
( v < Vpeak), on the other hand, increases as a function of in- 
creasing frequency. For example, the specific intensity mea- 
sured toward the Galactic polar region rapidly drops below 
3-5 MHz (where angular resolution is exceedingly poor), 
falling by an order of magnitude as the frequency is lowered 
to ^ 400 kHz, and by another orde r of magnitude by ^ 250 
kHz ( Brown 1973 ': 'Dulk et al. 2001). 

At relatively high frequencies, around 200-400 MHz, the 
spectrum steepens from spectral indices in the range 0.2-0.7 
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Table 4. Summary of low frequency ground-based radio observations" 
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"Data extracted from maps produced by Ellis & Mendillo 1987 (1.6 MHz), Ellis 1982 (2.1-16.5 MHz), Roger et al. 1999 (22 MHz), Alvarez 1997 (45 MHz), 
Landecker & Wielebinski 1970 (150 MHz), and Haslam et al. 1982 (408 MHz). 
''Symbols: u - frequency (MHz), R - angular resolution (degree). 

"^Brightness (in 10"'^ erg s"' cm"^ sr"' Hz"') measured toward various parts of the sky for which multi-frequency data exist: SP - the south Galactic pole 
(/,fo)=(0°,-90°); SPR- a south polar region (30°, -70°); GC - the Galactic center (/,^;)=(0° , 0°); NGC - the Galactic plane, near the Galactic center (20°, 0°); 
NAC - the Galactic plane, near the Galactic anti-center (240° , 0° ); and LI - a region of low brightness within loop I (320° , 50° ). 



to s ~ 0.8-0.9, with small deviations in the sky, the largest 
of which are associated with the Galactic spurs (loops I, III 
and IV, which exhibit spectra flatter in their center and steeper 
in their ridges, see Bridle 1967; Webster 1974). At very low 
frequencies, below 300 kHz, an excess radiation has been ob- 
served (Brown 1973; Manning & Dulk 2001), flattening the 
spectrum. Although the excess was previously accounted for 
by a hot, tenuous ISM phase (Fleishmann & Tokarev 1995), 
recent observations indicate that the excess radiation is lo- 
calized (at ecliptic coordinates led — 30° or led — 210°, see 
iManning & Dulk 2001), suggesting that the radiation origi- 
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Fig. 8. — Low frequency observations. Ground-based data (error bars, 
lines are aguide to the eye) are shown for five regions in the sky (see text 
and Tablegl: the Galactic center [GC, (Lb) = (0,0), solid line]; near the GC 
[NGC, (20°, 0), dotted line]; near the anti-center [NAC, (240°, 0), dashed 
line]; a south polar region [SPR, (30°, -70°), heavy line]; and a low bright- 
ness region near loop I [LI, (320°, 50°), dash-dotted line]. In spite of the 
large statistical and systematical errors (involving poor and frequency depen- 
dent angular resolutions, temperature scales, and zero-point temperatures), 
the frequency-dependent anisotropy pattern is apparent. The low resolution 
satellite data (see legend) approximately reflect the brightness at high Galac- 
tic latitudes, which dominate the sky at very low frequencies. 



nates from a local source. We shall not deal with low fre- 
quencies below 400 kHz in the following discussion. 

Anisotropy Pattern 

The LFRB is in general anisotropic, the radio maps re- 
vealing structure on various scales. Structure on hemispheric 
scales is evident, indicating that most of the radiation is Galac- 
tic, and thus reflecting the structure of the Galaxy and the po- 
sition of the solar system within it. The anisotropy pattern of 
the LFRB is highly frequency dependent. Figure |8l illustrates 
this anisotropy pattern, by comparing the spectra measured 
towards different directions in the sky. 

All experiments indicate that for frequencies v ^6 MHz, 
the background is strongest toward the Galactic plane, the 
brightest emission originating from the direction of the Galac- 
tic center itself. In the frequency range 600 kHz < v < 
2.5 MHz, on the other hand, the radiation is dominated by 
emission originating from the Galactic polar regions, ac- 
cording to most of the experiments sensitive to this en- 
ergy range: measurements by the RAE-1 satellite indicate 
that the radiation from the polar region is stronger than 
the radiation fr om the Galactic center and anti-cente r be- 
low - 1 MHz ( Clark eta l. 197(1 [Alexander et al. 197(1) : the 
IMP-6 experiment suggests with a high probability that the 
radiation between 130 kHz and 2.6 MHz is maximal to- 
wards the Galactic pol es and minimal towards the eclip- 
tic plane iBrown 1973h : and the WAVES experiment in- 
dicates that the radiation is dominated by emission from 
the Galactic polar regions between ~ 600 kHz and ^2.5 
MHz, is essentially isotropic at ^ 3.5 MHz, and is domi- 
nated by emissio n from the Galactic center above ^ 5 MHz 
(Manning & Dul k 20011) . On the other hand, it is interesting 
to note that observations carried out by the RAE-2 space- 
craft find no such phenomenon and indicate domination of 
the Galactic disk over the Galactic polar re gions at all fre- 
quencies above ~ 500 kHz (Novac o & Brown 1978h . Most 
illuminating are the sky maps obtained by very low fre- 
quency (1.6-5.5 MHz) ground-based telescopes positioned 
in Tasmania (Ellis 1982; Ellis & Mendillo 1987), which sug- 
gest a transition in the anisotropy pattern taking place between 
2.1 MHz and 3.7 MHz; whereas above 3.7 MHz the emission 
is stronger toward low latitudes (near the Galactic plane), the 
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pattern below 2.1 is reversed, stronger emission detected to- 
ward high latitudes (near the Galactic poles). This effect, al- 
though small, is significant and consistent with most space- 
based observations carried out at low frequencies and with all 
ground-based observations at high frequencies. 

4.2. Galactic Model 

Since the 1950s, it has been widely accepted that the 
LFRB is dominated by Galactic synchrotron radiation that 
is emitted as relativistic (cosmic-ray) electrons gyrate in 
the magnetic fields reposing the ISM. At higher fre- 
quencies, the synchrotron signal is overtaken by compet- 
ing radiative processes, namely the CMB (above 1 GHz), 
thermal bremsstrahlung (d ominant above several GHz, 
see e.g. iLawsonetal. 198% . and emi ssion from spinning 
dust (dominant above ~ 70 GHz, e.g. iTe gmark et al. 200(1 
iKnkbeiner 2003). Transition radiation, emitted when charged 
particles are accelerated in a medium with a varying refrac- 
tive index, dominates at very low frequencies (i^ <C 100 kHz, 
^leishmann & Tokar ev 19951) . An additional, extragalactic 
radio component must also be present, at least partly at- 
tribute d to the inte grated emission from discrete radio sources 
(e.g. Simo n 19771 and the references therein) and to inter- 
galactic shocks. 

In this section we construct a simple Galactic mod el, ba sed 
on important features of the observed LFRB (see ^4. II and 
Figure |8} and the present understanding of the Galactic dis- 
tributions of cosmic-rays, ionized gas and magnetic fields. 
Large observational uncertainties and poor resolution render 
all but the simplest qualitative model redundant, but the lat- 
ter is sufficient to show that the observed background may, 
in principle, be entirely Galactic. The good fit of our Galac- 
tic model to existing observations (within present uncertain- 
ties) shows that the introduction of an additional, extragalactic 
component is not necessary for explaining the observations. 

Construction 

We begin by considering an over-simplified model, in 
which the Galactic composition is assumed uniform along a 
given line of sight; more specifically, we assume that the emis- 
sivity j(i^) and the absorption coefficient aiv) are roughly 
constant throughout the Galaxy, along the line of sight. Sup- 
plementing the Galactic emission by an additional extragalac- 
tic component of specific intensity /i(j^), thus leads to a simple 
two-component model. The equation of radiative transfer may 
then be integrated to give the specific intensity as a function 
of the distance L transversed along the line of sight through 
the Galaxy: 



t(v) 



(31) 



where T(iy) = a(v)L is the optical depth. In the limits of very 
small or very large optical depths, this expression simplifies 
to: 

i(„ n ~ j i{v)L+IAv) foxTiv) < 1 ; 

iKV.L^) - <y j(,y)/Q,(jy) + /^.(i.)e-^M for t{v) > 1 . 

Note that such a model implies that the observed intensity is 
a monotonically rising function of L for frequencies where 
I.xiv) < i{v)/a{v) and is nearly isotropic in low frequencies, 
in contrast to the anisotropy pattern observed. In addition, 
uniform emission and absorption can not explain the observed 
sharp cutoff at frequencies v < 800 kHz (although for certain 



ISM parameters the cutof f may be attributed to the R azin- 
Tsytovich effect, see e.g. iFleishmann & Tokarev 19951 and 
the references therein). 

The observed low-frequency cutoff at low Galactic latitudes 
scales roughly as cx u'^ with 5 ~ -2 (see Figure char- 
acteristic of mixed emission and absorption at a high optical 
depth [see equation ( I32H . On the other hand, at high latitudes 
the cutoff is stronger, i < -2.7 for ly < 800 kHz (in contrast, 
for ly >20 MHz the spectrum is typically softer at high lati- 
tudes). Such a cutoff is naturally explained by the presence 
of a thin layer of relatively dense, strongly absorbing ionized 
gas near the Galactic plane. A double-component Galactic 
model, consisting of a thin, dense disk and a thick, dilute disk 
(as found by previous models which have focused on the high 
frequency angular data, see ^4.3t . may also naturally explain 
the observed frequency-dependent anisotropy pattern of the 
LFRB. Unlike a single-component Galactic model (with or 
without an additional isotropic, extragalactic component), in 
such a model the brightness need not be a monotonic func- 
tion of the distance L transversed through the Galaxy along 
the line of sight [for frequencies where /i(i') < j(i')/a(i')]. 

Finally, the spectrum measured at a given direction is often 
well fit at low (ly < few MHz) and at high (i/ > few 10 MHz) 
frequencies by equation i32\ . Hence, one may directly ex- 
tract from the data the average values of i{v)/a{v) and i{v)L 
towards various lines of sight [note that equations and 
(|32} should be slightly modified for a two-component Galac- 
tic model]. These values lead to a set of equations, relating the 
parameters of the Galactic model. Breaking the degeneracy of 
the model with some assumptions than yields the parameters 
in a straightforward fashion. 

Results 

Figure |5] illustrates how a simple Galactic model can ac- 
count for the observed frequenc y-de pendent anisotropy pat- 
tern of the LFRB, discussed in ^4.11 Since the examined re- 
gions on the sky are only representative and the uncertainties 
are large, we can only demonstrate a qualitative agreement 
with the data, and argue that the model captures the important 
observed features. 

The model consists of two uniform Galactic disks; a 
thin disk (denoted by subscript 1), containing most of the 
ionized gas, and a thick disk (denoted by subscript 2). 
The small number of sky regions examined and the un- 
certainty and low resolution of the data, imply that the 
model parameters are highly degenerate. Demanding that 
the 10 MHz synchrotron emissivity at the Galactocentric ra- 
dius of the solar system (/?, ~ 8 kpc) agrees with screening 
method measurements ( /V ~ 3 x lO"-'^ erg s"' cm"-' sr"' Hz~', 
see Fleishmann & Tokarev 1995, and the references therein), 
helps break the degeneracy, yielding the best fit parameters 
used in Figure|9] The thin disk is found to have a radius R\ ~ 
8.5 kpc, a half scale height z\iRs)/2 ~ 20 pc at the Galac- 
tocentric radius of the solar system, synchrotron emissivity 
ji(iy) ~ 6.6 X 10-4'V/10MHz)-°-25 gj.g ^-i ^,^-2 ^^^-i 

and an absorption coefficient ai(l MHz) ~ 1.5 x 10"^' cm"'. 
For the thick disk, we find R2 ~ 10 kpc, ziiRs)/'^ — 1 kpc, 
j2{v) ~ 2.4 x 10-^'^(i^/10MHz)-° ''5 erg s"' cm'^ sr"' Hz"' 
and a2(l MHz) ~ 1.4 x lO'^' cm"'. 

The dimensions of the Galactic disks found above are 
in qualitative agreement with concurrent Galactic models, 
based on various observations that measure emission and ab- 
sorption processes in the ISM, and propagation effects of 
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radiation p assin g throug h it (see th e discussion of spatial 
models in Phillips et al. 1981 at [Cordes & Lazio 20021 
iKeshet etal. 200 4a). However, the model parameters should 
be considered as order of magnitude estimates only, because 
they depend upon the modelled regions (chosen by the avail- 
ability of multi-frequency data) and on the choice of degen- 
eracy breaking. In addition, these values represent averages 
over large, ~ kpc distances, where large fluctuation exist. 

For example, the emissivity of the thick disk is found to be 
a factor of 3.5 higher than the emissivity of the thin disk for 
I' =10 MHz (the two become comparable for i' ~ 240 MHz). 
However, this reflects the sensitivity of the high latitude ob- 
servations modelled to the emissivity in the Galactocentric ra- 
dius of the sun , known to be enhanced with in the local arm 
dCaswell 19761 'Fleishma nn & Tokarev 1 995). and does not 
imply that the emissivity is overall larger in the thick disk. 
Similarly, in our simple model, the value of a is strongly 
constrained by measurements in low frequencies, where the 
optical depth is large and the radiation observed originates 
from the vicinity of the sun. Hence, the modelled value of a 
reflects the state of the ISM in the solar neighborhood, and 
is unsurprisingly found to be similar in both disks. The ab- 
sorption coefficient provides a measure of the thermodynam- 
ical state of the modelled ISM, roughly proportional to the 
value of the combination n^T'^l^ averaged over the differ- 
ent phases of the ISM along the line of sight, where T is the 
temperature and iig is the electron number density. For typi- 
cal values a(l MHz) ~ 1.5 x 10"^' c m"' and a local density 
Ue ~ 0.05 cm-3 (e.g. Revnolds T99ll) . we find {T) ~ 2500 K. 
These estimates for the average values of J{v), a{v) and T 
may not be extrapolated to large distances away from the so- 
lar system based on the present analysis. 

Our model is highly oversimplified, lacking for example 
spiral structure, radial gradients along the Galactic plane, and 



north-south Galactic asymmetry. It is important to note, how- 
ever, that within observational uncertainties, the model pro- 
vides a good fit to the data without the introduction of an 
additional, isotropic extragalactic component. Hence, the 
data does not directly imply the existence of a DERB. The 
frequency-dependent anisotropy pattern of the LFRB may be 
used in the future, not only to construct more realistic Galac- 
tic models, but also to measure the extragalactic radio back- 
ground or to place stringent limits on its magnitude. For 
example, improved data may allow one to study the bright- 
ness of the polar regions at low optical depths, and to mea- 
sure the deviation from an angular disk-like profile or from a 
pure power-law spectrum. With the present data uncertainties, 
however, we can only impose upper limits on the DERB flux 
which are of the same order as the measured Galactic signal. 

4.3. Previous models 

In the following we present the main LFRB models that 
have appeared in the literature. Although various models at 
different levels of sophistication have been proposed to ac- 
count for either spectral or angular features of the observed 
LFRB, no comprehensive model has thus far accounted for 
the combined spectral and angular data. In particular, we are 
not aware of any model that explains the frequency-dependent 
anisotropy pattern of the LFRB discussed above, although the 
basic concept of enhanced absorption at low latitudes has been 
acknowledged long ago (e.g. Ellis 1982). In general, models 
accounting for the spectral features observed (hereafter spec- 
tral models) yield little information about the structure of the 
Galaxy, but can in principle separate between a Galactic and 
an extragalactic component. Models accounting for the an- 
gular data (hereafter spatial models) can be used to study the 
structure of the Galaxy, but do not provide information re- 
garding the extragalactic component; the latter must therefore 
be separately specified and subtracted from the data. 
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Fig. 9. — Demonstration of a simple model which reproduces the qualita- 
tive features of the observed LFRB anisotropy pattern, different lines corre- 
sponding to different regions in the sky. The model consists of two Galactic 
disks (thin d isk a nd thick disk) and no extragalactic component. For model 
details, see i|4.2l Notations for ground-based observations (see legend) are 
the same as in Figure |8] but here the lines represent the predictions of the 
model. The minimal and maximal brightness measured by the IMP-6 space- 
craft (triangles) are for integrated large fields of view, and should thus be 
modelled by a frequency-dependent sky averaging scheme. 



Spectral Models 

Spectral models have been used to explain the radio spec- 
trum measured towards a given beam direction, without ac- 
counting for the anisotropy pattern. Such models are es- 
pecially appropriate for low resolution LFRB observation, 
such as satellite measurements, which provide mostly spec- 
tral information. Spectral models often make two simplify- 
ing assumption: (i) the radio signal is composed of Galactic 
foreground and an extragalactic background, represented by 
the (oversimplified) two-component model discussed in 14.21 
[see equation ( I31H : and (ii) both Galactic and extragalactic 
components are assumed to have intrinsic power-law spec- 
tra, jiv) oc v'^K and Ix{v) oc v~^', at low frequencies (say, 
V < 200 MHz). With these assumptions, one may use devia- 
tions of the observed intensity from a pure power-law, in order 
to fit the parameters of the model. 

In this method, Clark et al. (1970) analyzed the emis- 
sion from the 'north halo minimum', a region of minimal 
radio brightness around Galactic coordinates / = 150° and 
b = 50°, which features a slight excess of emission in frequen- 
cies 2-4 MHz. Attributing this excess to extragalactic emis- 
sion yields Sg ~ 0.4±0.05, s, ~ 0.8±0.1 and 7,(10 MHz) ~ 
(3.5 ± 1.3) X 10"'** erg s"' cm"^ sr"' Hz~'. However, in or- 
der to fit the data, the extragalactic background must cut off 
sharply below 3 MHz, presumably because of strong extra- 
galactic absorption. Cane (1979) has analyzed the emission 
from the Galactic poles, and reported that although the data 
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may be explained by a uniform Galactic disk, a better fit is ob- 
tained by including an extragalactic component. By assuming 
that s, ~ 0.8, she finds Sg ~ 0.51 ±0.02 and 7,(10 MHz) ~ 
(1.6 ±0.5) X 10"^^ erg s~' cm~^ sr"' Hz"'. The synchrotron 
emissivity has been estimated by attributing the high lati- 
tude radio brightness to the integrated emissivity through a 
uniform disk of some assumed half scale height z- Thus, 
Cane (1979) assumed z = 500 pc, and found 7(10 MHz) ~ 
(4.9 ±0.3) X 10-^''erg s"' cm"^ sr"' Hz"'. 

In addition to the shortc omin gs of the above oversimpli- 
fied model, as discussed in ^4.21 such spectral models may in 
general lead to an overestimated extragalactic low-frequency 
radio background. The reason for this is that the modelled 
Galactic emission alone provides an acceptable fit to the spec- 
tral data, such that the additional extragalactic component 
must be smaller than the Galactic component at all frequen- 
cies and in all directions. Hence, the reported extragalactic 
background was found to be on the order of the expected in- 
tegrated emission from discrete radio sources (Cane 1979), 
or some what stronger bu t intrinsically cutting off below ~ 
3 MHz JClark et al. 197(1) . The preceding discussion indi- 
cates, that without a significant spectral feature, that persists 
over various lines of sight, and considering the present obser- 
vational uncertainties and limited understanding of Galactic 
radio emission, spectral models can not reliably identify the 
extragalactic component. Spectral features observed at very 
low frequencies < 10 MHz), where the spectrum of Galac- 
tic emission strongly deviates from a power-law, may in prin- 
ciple be used to estimate the DERB flux, but (i) with present 
uncertainties such estimates can only be used as upper limits; 
and (ii) these limits are only applicable for very low frequen- 
cies, say < 30 MHz. Independent attempts to estimate the 
extragalactic emission by making assumptions regarding the 
uniformity of the Galactic spectral index Sg (e.g. Bridle 19671), 
are sensitive to this inaccurate assumption and may misinter- 
pret an isotropic Galactic component as being extragalactic, 
as pointed out by Lawson et al. (1987). 

Spatial Models 

More sophisticated m odels of Galactic radio emission 
iPhilUps et al. 19 8Ta: Beuerm ann et al. 1 985) have been con- 
structed by analyzing detailed all-sky maps at frequencies 
where absorption is negligible, su ch as the Bonn all- sky ra- 
dio continuum survey at 408 MHz (Haslam et al. 1982). Such 
three-dimensional models generally identify the distribution 
of Galactic cosmic-ray electrons and magnetic fields, as com- 
posed of a disk or a combination of disks, with spiral structure 
directly evident from steps in the radio temp erature profile 
observed along the Galactic plane jMills 19591) . Unfolding 
procedures, whereby the brightness in a given line of sight is 
related to the various spiral sections contributing to it, were 
used to calculate the locations and the relative emissivities of 
the Galactic spiral arms. While the negligible absorption con- 
siderably simplifies the analysis and allows one to reconstruct 
the structure of the Galaxy, it is difficult to distinguish be- 
tween the Galactic signal and the extragalactic background. 
Spatial models thus assume (and remove) some isotropic ex- 
tragalactic component. For example, a component of bright- 
ness temperature ~ 6 K is often uniformly subtracted from the 
408 MHz temperature map in order to account for extragalac- 
tic emission, assumed to be composed of similar contributions 
from the CMB and from discrete radio sources. 

As an example, Phillips et al. (1981a) have unfolded the 



emission originating from the Galactic plane, finding spiral 
structure at distances between 3.6 kpc and 16-20 kpc from 
the GC with at least three spiral arms of an arm-inter-arm 
emissivity ratio > 5, and placing the solar system between 
two major arms. The radial emissivity profile within a given 
spiral arm is degenerate in such models. Hence, a good fit 
to the data was found for various profiles, including an ex- 
ponential profile of scale length 3.9 kpc, provided that the 
Galaxy contains roughly equal magnitudes of regular (ori- 
ented along its spiral structure) and random magnetic fields. 
Extending the analysis to the entire sky (Phillips et al. 198 Itl) 
indicates that ^ 90% of the emission originates from a disk 
of HWHM thickness 0.46 -0.85 kpc at the Galactocenti'ic 
radius Rg of the sun, and ~ 10% originates from a much 
thicker disk, extending ^ 10 kpc from the Galactic plane. 
Both disks grow thicker at larger distances from the GC, re- 
sembling the thi ckening of the atomic hy drogen disk outside 
the solar radius iJackson & Kellman 1974) . A good fit to the 
data was obtained when the thinner disk scales as ~ exp(-5 + 
R/IO kpc). The implied emissivity at the solar neighborhood 
is ;(408 MHz) ~ (0.4- 1.9) x lO"'*" erg s"' cm'^ sr"' Hz"', 
suggesting the presence of regular and irregular magnetic 
fields of magnitude 3-6 /iG each. Extrapolating a typical 
emissivity ;(408 MHz) ~ 1.5 x 10""'" erg s"' cm"^ sr"' Hz"' 
to 10 MHz frequencies, by assuming an effective spectral in- 
dex s,„ = 0.6 in the range 10-408 MHz, gives ;(10 MHz) ~ 
1.4 X 10"^^ erg s"' cm'^ sr"' Hz"'. 

Clearly, such single-frequency studies are limited in their 
capability to reconstruct the structure of the Galaxy, as evi- 
dent from the degeneracies of the model (e.g. in the radial 
profile) and some of its non-physical features (e.g. the ex- 
aggerated scale height of the thick disk). By performing a 
detailed spatial analysis using such unfolding procedures, but 
for sky maps measured in more than one frequency, in partic- 
ular at low frequencies where absorption is not negligible, one 
may construct a far better model of Galactic emission and ab- 
sorption than currently available. The combined spectral and 
angular data will remove some of the degeneracies present in 
both spectral and spatial models, relax some of their simpli- 
fying assumptions, and possibly provide a direct estimate of 
the extragalactic component. Such an approach requires low 
frequency sky maps of better resolution than presently avail- 
able, and is beyond the scope of this paper. However, future 
very low {i/ < 10 MHz) frequency observations, for example 
by the ALFA mission, may enable such an analysis. 

5. DISCUSSION 

We have studied the radio signal produced by synchrotron 
emission from the strong intergalactic shocks associated with 
structure formation (see |2j. The analytical model of Wax- 
man & Loeb (2000) was generalized by adapting it for a 
ACDM universe, incorporating spectral features and shock 
asymmetry into the model, and calibrating the free parame- 
ters of the model, as summarized in Table |2l The halo param- 
eters face, fr, and fr were calibrated using a hydrodynamical 
cosmological simulation, by demanding that the model agrees 
with the simulation on the average baryon temperature, the av- 
erage mass consumption rate by strong shocks, and the typical 
size of regions where the thermal energy is injected by shocks. 
After calibrating the model with these essentially global fea- 
tures, it yields ra dio (see Figures |2l |3] and |5} and 7-ray (see 
iKeshet etal. 2003 ) signals which are in good agreement with 
the signals extracted independently from the simulation. Al- 
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though the parameter calibration scheme is inaccurate and the 
weak redshift dependence of the parameters needs yet to be 
determined by a more careful analysis, this agreement sug- 
gests that the calibration procedure is sensible. The local- 
ized nature of regions where most of the thermal energy is 
injected (e.g. at the intersections of X-ray cluster accretion 
sho cks with galaxy fi laments, channelling gas into clusters, 
see lKeshet et al. 200 3). enhances the synchrotron luminosity 
of a halo while leaving its inverse-Compton luminosity un- 
changed, because of the enhanced magnetic energy density. 
The radio luminosity-temperature relation according to the 
calibrated model was shown to be in qualitative agreement 
with observations of cluster radio halos (see Figure|4j. 

We have examined the observational consequences of the 
predicted radio signal, for present-day and for future radio 
telescopes (see Our main results are illustrated in Fig ures 
|6land0(for the calibrated model parameters discussed in ^2.31 
and summarized in Table|2j- Figure|6lshows the contribution 
of various signals to the radio sky, suggesting that emission 
from intergalactic shocks contributes up to a few tens of per- 
cent of the extragalactic radio background below 500 MHz 
(see also Figure depicts the angular power spectrum 
of various signals on an angular scale of ~ 0°.5, along with 
the sensitivities and the angular resolutions of the telescopes. 
The figure demonstrates that the designs of next generation 
radio telescopes such as the LOFAR and the SKA, are more 
than sufficient for a detection of the angular fluctuations in- 
troduced by intergalactic shocks. In fact, even present-day 
high-resolution radio telescopes (such as the VLA in its maxi- 
mum sensitivity configuration), are potentially sensitive to the 
predicted signal. 

Foreground fluctuations in the synchrotron emission of our 
Galaxy limit a positive detection of the signal to sub-degree 
scales and to high Galactic latitudes. Confusion with discrete 
radio sources requires that the brightest sources be modelled 
and removed. With a feasible point source cut (100 mJy at 
1.4 GHz or equivalently 3 Jy at 150 MHz), the intergalactic 
signal dominates over all foreground and background signals, 
at angular scales 10' < < 1° a nd fr equencies < 10 GHz 
(assuming ^g^^ ^ 3 x 10"^, see ^2.41 and the discussion be- 
low). For the above cut, the signal is of the same magnitude 
as the integrated emission from discrete sources at ~ 1' scales. 
Attempts to detect the signal at < 1' scales are confusion lim- 
ited by discrete sources, demanding increasingly more ambi- 
tious point source removal schemes. In addition to the labor 
involved, noise is always introduced when removing discrete 
sources, because of source model uncertainties. The spectrum 
of the signal suggests that it is most pronounced at frequen- 
cies around ^100 MHz, planned to be covered by both the 
LOFAR and the SKA. 

The calculated level of angular fluctuations in the radio 
emission from intergalactic shocks is sensitive to uncertain- 
ties in the calibration of the model parameters. The logarith- 
mic contribution to the variance (shown in Figure scales 
roughly according to 

51,(1- ^,,^B,faccjTjr) = Ub^^ Sllifrt, 1,1,1,1,1). 

Jr 

(33) 

The radio signal calculated from the ACDM simulation 
(jKeshet et al. 2004b), however, depends only on the en- 
ergy fractions and ^b, through the combination £,e^B- 
The results of the simulation, as shown in Figures |2] |3l 
and I5I7I suggest that Sli was in fact slightly underesti- 



mated by our parameter calibration scheme. This implies, 
for example, that intergalactic shocks introduce inten- 
sity fluctuations of magnitude Sli > 8 x 10~^''\Ce / 5 x 
10"'*)(i^/100MHz)"' erg s"' cm"^ sr"' Hz"' on multipoles 
400 < I < 2000, corresponding to temperature fluctuations 
ST, > 260(^,^b/5 X 10-4)(i//100 MHz)-^ K. 

Our model implies that emission from intergalactic shocks 
on angular scales ~ 0° .5 dominates the sky at < 500 MHz 
frequencies if ^^^b ^ 10""*, and is even dominant at frequen- 
cies as high as 10 GHz if S^^S^b ^ 3 x 10~^. The simulation 
suggests that the signal is stronger than predicted by the (cal- 
ibrated) model, in particular on small angular scales. There- 
fore, according to the simulation, emission from intergalactic 
shocks will dominate the sky at < 0°.5 scales and < 500 MHz 
frequencies, even if ^^^b — 5 x 10"^. We have used observa- 
tions of SNR shocks and of magnetic fields in the halos of 
galaxy clusters, in order to show that for strong intergalac- 
tic shocks S^eS^B — 5 X 10""*, unlikely to be smaller than this 
value by more than a factor of ~ 4. For this purpose, in ^2.41 
we have used dimensional analysis arguments to show that 
the physics of strong intergalactic shocks is essentially iden- 
tical to the physics of strong SNR shocks, both of compa- 
rable velocities v ~ 10^ km s"', provided that an appropriate 
re-scaling of time is carried out (see also Keshet et al. 2003). 
We have assumed that strong shocks accelerate electrons to a 
power law distribution of index p = 2 (equal energy per loga- 
rithmic interval of electron energy). Such a distribution is in- 
ferred from SNR observations and agrees with linear models 
for diffusive shock acceleration, although non-linear models 
suggest deviations from a pure power-law. 

As noted by Waxman & Loeb (2000), emission from in- 
tergalactic shocks may have akeady been detected by CMB 
anisotropy studies at frequencies < 10 GHz and angular 
scales smaller than a degree. In particular, telescopes operat- 
ing at frequencies of a few GHz such as ACSA, Ryle and the 
VLA, may have detected the signal at arcminute scales. High 
resolution ~ 1 GHz radio surveys may have also detected the 
signal by now at an identifiable level. Such surveys are gener- 
ally available at low Galactic latitudes, where contamination 
by Galactic foreground fluctuations is severe. Nevertheless, 
analysis of 'quiet regions' observed off the Galactic plane, in- 
cluding a careful removal of discrete sources, may yield the 
desired signal. It will probably be easiest to identify the sig- 
nal at 1'- 10' scales, by modelling discrete sources and cross 
correlating the maps with known tracers of large-scale struc- 
ture. Detection of emission from intergalactic shocks is not 
unrealistic even in case the signal has been mildly overesti- 
mated. A signal dli lower than calculated in ^ by a factor 
of a few, may still be identified at ^ 1'- 10' scales, if faint 
point sources are modelled and removed. If the signal was 
overestimated by ~ 2 orders of magnitude, it may still be 
detectable by next generation radio telescopes, by means of 
cross-correlation with known tracers of large-scale structure 
such as galaxy counts, or, in the future, with 7-ray emission 
from intergalactic shocks. 

As mentioned in ^ future detection of radio emission from 
intergalactic shocks will have important implications on our 
understanding of cosmology and astrophysics. Detection of 
the signal will provide the first identification of intergalactic 
shocks, revealing the underlying cosmological flows and pro- 
viding a test for structure formation models. The signal, in 
particular when combined with 7-ray detection, will provide 
a measure of the unknown magnetic fields in the intergalactic 
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medium. Although non-trivial for interpretation, such a mea- 
sure of the magnetic field may provide insight into the un- 
known processes leading to IGM magnetization. The signal 
may also confirm the existence of the undetected warm-hot 
intergalactic medium, and provide a tracer for its distribution. 

Synchrotron emission from intergalactic shocks is corre- 
lated with the large-scale structure of the low-redshift (z < 1) 
universe, tracing young galaxy clusters and filaments. The 
signal could thus account for some observed features of the 
radio signature of galaxy clusters, namely radio halos and 
radio relics. Extended accretion shocks could contribute to 
radio halos (e.g. to the large rad io halo of the Coma super- 
cluster, see lThierbach et al. 2003L and the references therein), 
whereas localized shocks (e.g. where galaxy filaments chan- 
nel large amounts of gas into the cluster regions) may be re- 
sponsible for some radio relics observed at the outskirts of 
clusters [such as th e prototype relic found in the Coma super- 
cluster, 1253H-275 jEnBlin et al. 199 8). and the large-s cale ra- 
dio arcs ob served in clus ters A3667 (Rott2erin2 et al. 1997h 
and A3376 (Bagchi 2002*)], and perhaps also for some of the 
anomalous features observed in radio halos (e.g. in the unre- 
laxed clusters described by Govoni et al. 2004). 

Radio emission from intergalactic shocks is an important 
source of contamination for other radio signals, such as the 
low frequency CMB, Galactic synchrotron fluctuations on 
sub-degree scales, and competing extragalactic radio signals 
such as bremsstrahlung from Lya clouds, and 21 cm tomogra- 
phy (Zaldarriaga et al. 2003; Morales & Hewitt 2003; Loeb & 
Zaldarriaga 2003; Gnedin & Shaver 2003). For example, the 
integrated emission from intergalactic shocks introduces dis- 
tortions in the low-frequency average brightness temperature, 
at the level of 5T ~ 1.6(i^/3 GHz)"-' mK. Efforts to measure 
distortions in the CMB spectrum at low frequencies, e.g. in 
the next mission of the Absolute Radiometer for Cosmology, 
Astrophysics and Diffuse Emission (ARCADE**), may thus be 
sensitive to the signal at the low frequency bands (^ 3 GHz, 
see lFixen et al. 2 004: Ko gut et al. 20 04). We have shown that 
radio emission from intergalactic shocks constitutes a sig- 
nificant fraction of the extragalactic radio background below 
500 MHz. This component should thus be taken into account 
when evaluating the propagation of ultra-high energy photons 
with energies above lO'** eV, because at these energies the 
effect of the radio background on the transparency of the uni- 
verse is stronger than the effect of the CMB, and may thus 
be important in models for the origin of ultra-high energy 
cosmic -rays. 

We have analyzed observations of the diffuse low frequency 
radio background below 500 MHz (see and highlighted 
the frequency-dependent anisotropy pattern observed in fre- 
quencies 1 MHz "^v <\Q MHz (see Figure|8}. We presented 
a simple Galactic model, consisting of a thin disk, contain- 
ing most of the ionized gas, and a thick disk. This model 

' See|http://arcade.gsfc.nasa.gov/arcade/mdex.htnil| 



was shown to provide a good fit to the data (see Figure |9jl, 
considering the observational uncertainties. The model thus 
enabled us to (i) assess the feasibility of directly inferring 
the flux of the diffuse extragalactic radio background from 
presently available observations; and (ii) demonstrate how an 
elaborate Galactic model may be constructed in a straightfor- 
ward fashion, using the anisotropy pattern observed in very 
low frequencies. 

The agreement of the Galactic model with observations, 
and the large uncertainties and poor resolution of present low 
frequency observations, preclude a direct identification of the 
diffuse extragalactic radio background. At best, an upper limit 
on the extragalactic component may be imposed, at very low 
(< 30 MHz) frequencies where free-free absorption in the 
Galaxy becomes important. Such a flux constraint, a factor 
of a few lower than the Galactic foreground, is itself higher 
than the calculated signal from discrete sources and from in- 
tergalactic shocks by a factor of a few. This implies that (i) the 
true extragalactic background in ;/ ~ 10 MHz frequencies is 
probably smaller than the Galactic foreground by a factor of a 
few; and (ii) emission from intergalactic shocks constitutes a 
significant fraction of the total extragalactic signal, measuring 
up to a few tens of percent. 

We have reviewed previous models of the low frequency 
radio background, and demonstrated that previous claims for 
direct identification of the diffuse extragalactic component are 
based on insignificant, probably Galactic, spectral features. 
We pointed out that previous models fail to account for the 
combined spectral and angular data at very low (1 MHz <v < 
10 MHz) frequencies, where absorption is non-negligible and 
the anisotropy pattern of the sky is highly frequency depen- 
dent. Modelling the combined data can be used to disentan- 
gle the distributions of Galactic cosmic-rays, ionized gas and 
magnetic fields, and lead to an elaborate three-dimensional 
model of the Galaxy. Our model is too simplistic, and the 
data uncertainties too large, to provide reliable estimates of 
the physical properties of our Galaxy, other than confirming 
its double-disk structure and the presence of a ~ 1 kpc thick 
disk of cosmic ray electrons and magnetic fields, and impos- 
ing constraints on free-free absorption in the solar neighbor- 
hood. Future high resolution observations at very low fre- 
quencies, produced for example by the ALFA space mission, 
could thus provide valuable information regarding the struc- 
ture and the composition of the Milky Way. 
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